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Preface 


The  potential  usefulness  of  the  impact  damper  for 
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dynamics  and  mechanics  and  as  an  opportunity  to  acquire 
laboratory  experience,  but  also  as  an  opportunity  to  do 
basic  research  in  a  promising  field  of  vibration  control 

that  is  not  common  knowledge. 
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laboratories  or  journals  published  in  Europe,  Russia,  and 
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Abstract 

I 

~*The  usefulness  of  the  impact  damper  in  eliminating 
vibrations  is  studied  analytically  and  experimentally. 
Laboratory  models  of  vibrating  systems  are  constructed  to 
evaluate  the  performance  of  the  impact  damper  in  reducing 
or  eliminating  forced  and  free  vibrations.  A  computer  simu 
lation  of  a  single  degree-of-f reedom  primary  system  in  free 
vibration  employing  an  impact  damper  is  constructed  for  the 
same  purpose.  Laboratory  free  vibration  results  are  com¬ 
pared  to  the  computer  simulation  in  order  to  judge  its 
accuracy . 

The  computer  simulation  is  employed  to  determine 
the  impact  damper's  performance  in  free  vibration  as  the 
system's  parameters  are  varied.  Two  significant  measures 
of  the  damper's  effectiveness  are  obtained  as  approximate 
functions  of  the  system's  parameters. 

Observations  regarding  reduction  in  amplitude  and 
steady  state  motion  were  made  for  the  impact  damper  in 


forced  vibration. 


THE  USEFULNESS  OF  IMPACT  DAMPERS  FOR 


SPACE  APPLICATIONS 


I.  Introduction 


An  impact  damper,  also  referred  to  in  the  litera¬ 
ture  as  a  rattle  damper  or  an  acceleration  damper,  is  a 
simple,  passive  damping  device.  It  operates  by  allowing  a 
vibrating  primary  mass  to  go  through  a  series  of  collisions 
With  a  secondary  mass  carried  in  or  on  the  primary  mass. 
Figure  1  shows  one  of  the  simplest  models  of  an  impact 
damper,  with  a  primary  mass  M  free  to  travel  in  one  dimen¬ 
sion  only,  acted  upon  by  a  forcing  function  F(t) ,  a  secon¬ 
dary  mass  m,  a  spring  of  stiffness  k,  and  a  dashpot  with  a 
damping  constant  c. 

In  the  simplest  case,  the  motion  of  the  secondary, 
or  impacting,  mass  m  is  assumed  to  be  a  result  of  colli¬ 
sions  with  the  primary  mass  alone,  so  the  impacting  mass 
has  a  constant  velocity  between  impacts.  If  F(t)  is  sinu¬ 
soidal,  then  the  momentum  exchange  and  the  energy  dissipa¬ 
tion  resulting  from  the  impacts  usually  results  in  a 
decrease  in  the  amplitude  of  motion  of  the  primary  mass. 

If  the  primary  mass  is  in  free  vibration  (F(t)  =  0),  then 
the  impacts  cause  a  more  rapid  decay  in  the  amplitude  of  the 


motion  of  the  primary  mass  through  energy  dissipation  and 
momentum  exchange. 

This  simple  damper  could  be  of  practical  use  for 
space  applications  in  eliminating  the  unwanted  vibrations 
of  antennas,  telescopes,  or  any  other  flexible  structure 
which  tends  to -oscillate  about  its  intended  orientation. 

In  the  near  vacuum  of  space,  external  damping  forces  are 
essentially  zero.  Thus,  such  structures  must  have  internal 
damping  designed  into  them.  If  the  impact  damper  provides 
sufficient,  reliable  damping  without  adding  prohibitive 
mass  to  the  total  payload,  it  could  be  a  solution  to  some 
oscillation  problems. 

An  extensive  literature  search  turned  up  much  work 
on  the  effectiveness  of  the  impact  damper  in  reducing  forced 
vibrations.  Contradictory  conclusions  were  identified. 

Paget  (Ref  1)  probably  did  the  earliest  writing  on  the 
impact  damper,  but  the  first  serious  analytical  work 
appears  to  be  that  of  Lieber  and  Jenson  (Ref  2) .  In  their 
paper,  work  and  energy  considerations  were  used  to  solve 
for  the  one  degree  of  freedom  motion  of  a  primary  mass 
undergoing  perfectly  inelastic  collisions  with  a  secondary 
mass.  These  results  were  used  to  calculate  a  damping  fac¬ 
tor  which  was  verified  experimentally  through  comparison 
with  the  damping  observed  in  the  free  vibration  of  a  beam 
with  an  impact  damper  attached.  Their  solution  predicted 
the  impact  damper  would  be  most  efficient  (do  the  most  work 


per  cycle)  if  two  impacts  per  cycle  occurred  with  impacts 
equally  spaced  in  time. 

Grubin  and  Lieber  (Ref  3)  gave  a  more  straightfor¬ 
ward  solution  of  the  motion  of  the  system  for  collisions 
ranging  from  perfectly  inelastic  to  perfectly  elastic.  In 
Reference  (4) ,  it  is  shown  that  solutions  are  possible 
when  stable  and  symmetric  motion  is  assumed;  i.e.,  that 
two  impacts  occur  at  equal  time  intervals  during  the  cycle. 
This  is  referred  to  as  symmetric,  two  impact  per  cycle 
motion.  Such  motion  has  often  been  assumed,  and  in  Refer¬ 
ence  (5)  was  reported  to  occur  when  an  impact  damper  was  ' 
attached  to  a  cantilever  beam  in  forced  vibration.  Grubin 
and  Lieber  (Ref  6)  went  on  to  do  a  stability  analysis  on 
this  symmetric,  two  impact  per  cycle  motion.  Lieber  and 
Duffy  (Ref  7)  modeled  a  cantilever  beam  with  an  impact 
damper  as  a  system  composed  of  four  lumped  masses  and  used 
an  electric  analog  model  of  this  system  to  study  the  effects 
of  parametric  changes  on  the  dampers'  performance. 

Feygin  (Ref  8)  solved  and  did  a  stability  analy¬ 
sis  for  the  motion  of  an  impact  damper  similar  to  that 
shown  in  Figure  1,  but  with  the  motion  of  the  impacting 
mass  between  impacts  subjected  to  dry  friction.  Masri 
(Ref  9)  started  with  the  assumption  of  symmetric,  two 
impact  per  cycle  operation  and  solved  for  the  motion  of 
the  system  under  sinusoidal  excitation.  He  also  did  a 
stability  analysis  to  show  this  motion  did  exist  for  a 


wide  range  of  system  parameters,  and  verified  his  results 
experimentally  and  with  both  digital  and  analog  computer 
simulations . 

Masri  (Ref  10)  solved  for  the  motion  of  a  forced 
system  with  any  number  of  impacts  per  cycle.  Sadek  (Ref  11) 
assumed  two  impacts  per  cycle  and  used  a  Fourier  series 
representation  of  the  impacting  forces  to  come  to  the  con¬ 
clusion  that,  in  general,  symmetric  impacts  do  not  occur, 
especially  for  system  parameters  leading  to  the  maximum 
reduction  in  amplitude  of  the  primary  system.  He  used  a 
laboratory  model  to  verify  his  results. 

Sadek  and  Mills  (Ref  12)  solved  for  the  motion  of 
the  system  in  forced  vibration  with  the  impact  damper 
affected  by  gravity,  while  Sadek  and  Williams  (Ref  13) 
then  provided  a  stability  analysis  on  these  results.  Sadek 
and  Thomas  (Ref  14)  solved  for  the  motion  of  a  system  in 
forced  vibration  and  with  the  secondary  mass  attached  to  a 
spring  and  influenced  by  gravity. 

Masri  and  Sadek  have  both  published  several  papers 
on  impact  dampers  with  carefully  solved  equations  of  motion 
and  stability  analyses.  The  only  significant  difference  in 
their  approach  is  that  Masri,  and  many  other  authors, 
modeled  the  impacts  as  being  of  infinitely  short  duration, 
thereby  causing  a  discontinuous  change  in  the  velocities, 
but  not  the  positions,  of  the  two  masses.  Sadek  uses  a 
Fourier  series  representation  of  the  impact  force  and 
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treats  impacts  as  being  of  short,  but  finite  duration. 
Masri's  (Ref  9)  solution  for  the  sinusoidally  forced  system 
shows  symmetric,  two  impact  per  cycle  motion  to  be  pos¬ 
sible  for  a  wide  range  of  system  parameters,  and  that  this 
kind  of  system  gives  the  maximum  reduction  in  the  primary 
system's  amplitude.  Sadek's  (Ref  11)  solution  is  for  a 
system  at  a  specified  r?  Lo  of  secondary  mass  to  primary 
mass  and  at  a  specified  forcing  frequency.  This  solution 
gives  only  one  value  for  the  gap  in  which  the  secondary 
mass  travels  which  will  give  symmetry,  and  this  gap  does 
not  give  maximum  reduction  in  amplitude.  Masri  (Ref  10) 
also  predicts  that  the  amplitude  decreases  with  an  increase 
in  the  ratio  of  the  secondary  mass  to  the  primary  mass, 
while  Sadek  (Ref  11)  says  that  there  is  an  optimum  value 
for  this  ratio,  and  that  increasing  it  beyond  this  point 
increases  the  system's  amplitude. 

Roy,  Rocke,  and  Foster  (Ref  15)  did  an  analytical 
and  experimental  study  of  the  impact  damper  in  the  center 
of  a  beam  in  benalng  vibration,  using  both  a  simply  sup¬ 
ported  beam  and  a  beam  with  both  ends  clamped.  They  used 
both  a  closed  form  solution  for  the  motion  of  the  beam 
between  impacts  and  a  discrete  mass  model  of  the  system  to 
do  numerical  calculations  of  the  motion  of  the  beam. 

These  numerical  solutions  were  verified  with  experimental 
results.  All  previous  researchers  did  the  analytical  work 
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assuming  a  rigid  bodied,  single  degree  of  freedom  primary 
system. 

Dokainish  and  Elmaraghy  (Ref  16)  did  a  computer 
simulation  of  an  impact  damper  and  produced  a  series  of 
curves  from  which  damper  performance  can  be  predicted  for  a 
given  set  of  system  parameters.  Yamada  (Ref  17)  solved  for 
the  motion  of  a  sinusoidally  excited  impact  damper  similar 
to  Figure  1,  but  with  a  piecewise  linear  spring.  Other 
solutions  to  the  forced  motion  of  different  impact  dampers 
can  be  found  in  References  (18),  (19),  (20),  (21),  and  (22). 

Yasuda  and  Toyoda  (Ref  23) ,  considered  the  useful¬ 
ness  of  an  impact  damper  in  reducing  the  free  vibration  of 
a  lightly  damped  system.  They  used  experimental  results 
to  obtain  parametric  relations  which  could  be  used  to 
solve  for  the  damping. 

The  purpose  of  this  thesis  is  to  examine  the  use¬ 
fulness  of  the  impact  damper  in  reducing  both  forced  and 
free  vibrations.  A  laboratory  model  and  a  computer  model 
of  a  freely  vibrating  system  with  an  impact  damper  were  con 
structed.  These  models  were  used  to  determine  damper  per¬ 
formance  as  coefficient  of  restitution  e,  mass  of  impacting 
mass  m,  distance  between  impacting  surfaces  d,  and  other 
parameters  were  varied.  The  forced  vibration  case  was 
examined  using  a  laboratory  model  consisting  of  an  upright 
flexible  beam  with  an  impact  damper  on  top  and  a  sinusoidal 
angular  displacement  applied  to  the  bottom. 


II .  Analytical  Studies 


The  motion  of  the  free  vibration  and  forced  vibra¬ 
tion  impact  damper  is  considered  in  this  chapter.  The 
analysis  is  for  the  types  of  impact  dampers  depicted  in 
Figures  2  and  3,  which  are  the  types  used  in  the  laboratory 
studies.  The  motion  of  the  free  vibration  impact  damper  of 
Figure  2  can  be  followed  analytically  through  any  number  of 
impacts.  The  solution  to  the  motion  of  the  forced  vibra¬ 
tion  impact  damper  of  Figure  3  is  not  completely  described 
in  this  thesis,  but  some  useful  information  is  obtained. 

The  equations  given  in  this  chapter  are  derived  in  detail 
in  Appendix  A. 

Free  Vibration  Impact  Damper 

The  primary  system  of  the  free  vibration  impact 
damper  of  Figure  2  is  the  damper  assembly,  which  provides 
the  impacting  surfaces,  and  the  beam.  The  impacting  mass 
is  not  considered  part  of  the  primary  system.  The  angular 
displacement  9  of  the  primary  system  can  be  described  as  a 
rotation  about  an  axis  perpendicular  to  the  plane  of 
Figure  2  and  containing  point  0 ,  called  the  rotation  point. 
Between  impacts,  the  primary  system  is  acted  upon  by  the 
flex  plate,  gravity,  and  viscous  damping.  The  equation 
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of  motion  of  the  primary  system  between  impacts  is,  for 
small  values  of  0 : 


10  +  c0  + 


El 


Wmr_  + 
T  Om 


0=0 


(1) 


where : 

I  =  moment  of  inertia  of  the  primary  system  about  0, 
with  units  of  mass  *  length^; 

c  =  viscous  damping  constant,  with  units  of  force  • 
length  •  time; 

WT  =  total  weight  of  the  primary  system; 

Tq  =  distance  from  0  to  the  center  of  mass  of  the 
L  primary  system; 

E  =  modulus  of  elasticity  of  the  flex  plate 
material; 

I*  =  moment  of  inertia  of  the  cross-sectional  area 
of  the  flex  plate,  with  units  of  length4;  and 

ds  =  length  of  the  flex  plate. 

This  is  more  conveniently  used  in  the  form: 


10  +  c0  +  k0  =  0 


(2) 


where : 


k 


Wmrn 
T  0m 


and  k  has  units  of  force  •  length.  The  motion  of  the  pri¬ 
mary  system  between  impacts  is  described  by  Equation  (1) 
if  the  rotation  point  0  is  stationary.  This  is  shown  to 
be  approximately  true  in  Appendix  A,  where  point  0  is  shown 


to  be  located  where  the  center  of  the  flex  plate  is  when 
it  is  undeformed. 

Solving  Equation  (2)  for  0  between  impacts  i  and 
i+1  gives: 


0  (At)  =  e 


<-ic  At) 


0(t  (  +  ))  0(t  (  +  ))c. 

- q -  +  2iq  1  Sin  qAt 


+  0 ( t^ ^ )  cos  qAt 


(3) 


where  At  is  the  time  since  impact  i,  0(t^  +  ^)  and  0(t^  +  M 
are  the  angular  position  and  angular  velocity  after  impact 
i,  and 
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If  there  is  no  viscous  damping  (c=0)  then: 


0  (+) 

0  (At)  =  -  sin  qAt  +  0(t.  )  cos  qAt 

q  l  ^ 


(4) 


and 


-Vi" 


The  equation  of  motion  for  the  suspended  impacting 
mass  between  collisions  is: 


2" 

mL  a  +  c  a  +  mgL  a  =  0 
m  a  m 


(5) 


where: 


m  =  mass  of  impacting  mass; 


L  =  distance  of  the  impacting  mass  from  its  rota- 
m  tion  point; 

a  =  angular  displacement  of  the  strings  suspending 
the  impacting  mass  from  the  vertical; 

ca  =  viscous  damping  constant,  with  units  of  force  • 
u  length  •  time;  and 

g  =  acceleration  of  gravity. 


For  small  values  of  a,  this  motion  is  almost  entirely  in 
the  horizontal  direction  or  x  direction.  Therefore,  Equa¬ 
tion  (3)  can  be  approximated  as: 


mx  +  cx  +  kx  =  0 
m  m 


(6) 


where : 


x  =  La; 
m 

2 

c  =  c  /L  ;  and 
m  am 


k  =  mg/L  . 
m  m 


c  has  units  of  force  •  time/length  and  k  has  units  of 
m  m 

force/length . 

Solving  Equation  (6)  for  x  between  impacts  i  and 
i+1  gives: 


x(At)  =  e 


«-25At' 


x{t.  ^  )  x(t.  ^  )c  , 

1  +  - i-rr - E)  sin  qmAt 


‘m 


2  mq 


m 


‘m 


+  x (t . )  cos  q  At 
x  m 


(7) 
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where  x(t^'  ')  and  x(t^  ')  are  the  position  and  velocity 
immediately  after  impact  i,  and: 


If  there  is  no  viscous  damping,  then: 
x(t.(+)) 

x(At)  =  - — -  sin  q  At  +  x(t.^)  cos  q  At  (8) 

q„  m  x  m 

m 

and: 

q 

^m  »  m 

Finally,  if  there  is  no  viscous  damping  or  gravity  effects, 
the  impacting  mass  has  constant  velocity  between  impacts, 
so: 

x ( At)  =  x(t±(+))  +  x(ti(+))At  (9) 


The  initial  conditions  used  in  the  solutions  to 

Equations  (2)  and  (6)  for  0  and  x  during  the  motion  between 

•  • 

impact  i  and  impact  i+1  are  9,  0,  x,  and  x  evaluated  at 

time  t^+^  (immediately  after  impact  i)  .  If  impacts  are 

assumed  to  be  of  infinitely  short  duration  (during  which 

•  • 

0  and  x  remain  unchanged,  and  0  and  x  change  discontinuously) 

then  the  initial  conditions  can  be  obtained  in  terms  of 
•  * 

8,  0,  x,  and  x  evaluated  at  time  t. 


(immediately  berore 


impact  i) .  The  positions  are  given  by: 


e(ti(+))  =  0(ti(  }) 

x(ti(  +  ))  =  xtt^"1) 

Angular  momentum  is  conserved  across  the  impact,  so: 

iett^"*)  +  mx(ti("))L  =  I0(ti(  +  ))  +  mx(ti(+))L 

Also,  the  velocities  across  an  impact  are  related  by: 

e[x(ti(_))  -  0(ti(_))L]  =  0(ti(+))L  -  x{ti(+)) 

where  e  is  the  coefficient  of  restitution.  These  two 
relations  give: 


(10) 

(11) 


x(t.(+))  =  — [0(t.  M  )I(l+e) 
1  I+mL 


+  x(ti(“))  (mL  -  ^)] 


(12) 


6(t.<+>) 


1 

I+mL2 


[0(ti(_)) (I-mL2e) 

+  x(t^~^  )mL(I+e)  ] 


(13) 


Using  Equations  (10)  to  (13),  the  position  and 
velocity  of  both  the  beam  and  impacting  mass  are  obtained 
at  time  t^  +  ^  (immediately  after  the  ith  impact)  in  terms 


L5 


of  the  positions  and  velocities  at  time  t^  (immediately 
before  the  ith  impact) .  If  the  overall  system  is  started 
with  known  initial  conditions  at  time  tg=0,  and  if  the  times 
of  impacts  t^,  t^...  are  known,  the  exact  solutions 

up  to  time  t^  ^  can  be  solved  in  terms  of  the  initial  con- 

(. f ) 

ditions  at  tg..  The  initial  conditions  at  t^  can  then 
be  solved  in  terms  of  the  final  conditions  at  t^  ^  ;  these 
initial  conditions  can  be  used  to  solve  for  the  exact  solu¬ 
tions  from  t^  +  ^  to  tj^  Initial  conditions  and  exact 

solutions  can  then  be  obtained  from  time  t2^  to  t.^ 
and  this  process  can  be  continued  for  as  many  inpacts  as 
is  desired. 

This  process  assumes  that  all  impact  times  t^  are 

known.  Actually,  impact  time  t i  must  be  found  by  itera- 

(+) 

tion,  using  the  known  solutions  after  time  t^_^  and  the 
requirement  that  the  suspended  mass  must  remain  between  the 
two  stops.  The  time  t^  is  then  defined  as  the  time  when 
the  impacting  mass  first  comes  in  contact  with  either  stop 
after  time  tj+|. 

This  solution  technique  was  used  in  the  two  com¬ 
puter  programs  of  Appendix  C.  The  first  program,  IDEAL, 
uses  Equations  (4)  and  (9)  to  trace  the  motion  of  the  pri¬ 
mary  system  and  impacting  mass  in  an  environment  with  no 
viscous  damping  or  gravity.  The  second  program,  LABSIM, 
uses  equations  (3)  and  (7)  to  include  viscous  damping  and 
gravity  effects  on  the  system.  Included  in  Appendix  C 


are  the  results  obtained  by  running  programs  IDEAL  and 
LABSIM. 

Program  IDEAL  was  run  for  a  variety  of  values  for  M, 
L,  e,  q,  and  d,  where  d  is  the  distance  between  the  impact¬ 
ing  surfaces.  The  results  of  one  of  these  runs  is  given 
in  Figure  4.  Whenever  m,  L ,  and  q  are  greater  than  zero, 
e  is  positive  and  less  than  one,  and  d  is  positive  and  less 

than  2  0„L  (d  <  2  6  ,  L  is  necessary  for  impacts  to  occur) , 
max  max 

the  impact  damper  will  initially  bring  about  a  rapid  reduc¬ 
tion  in  amplitude.  After  the  maximum  amplitude  attained 
by  the  primary  system  during  each  cycle  declines  to  a 
certain  value,  the  impact  damper  becomes  much  less  effec¬ 
tive,  with  a  much  slower  reduction  in  amplitude  of  the  pri¬ 
mary  system.  In  the  region  where  the  damper  is  effective, 
the  maximum  amplitude  attained  is  observed  to  decrease 
approximately  linearly  with  time,  and  a  linear  function 
is  fit  to  the  peaks  using  a  least  squares  method.  The  de¬ 
cline  in  maximum  amplitude  with  time  is  denoted  by  s,  and  is 
a  measure  of  the  impact  damper's  performance.  The  maximum 
amplitude  attained  when  the  impact  damper  becomes  almost 

ineffective  is  denoted  0  ,  and  is  also  a  measure  of 

mr 

damper  performance. 

It  is  important  to  determine  if  s  and  0^  are  depen¬ 
dent  upon  parameters  other  than  M,  L,  e,  q,  and  d.  All  of 
the  computer  simulations  of  the  impact  damper  began  with 
the  primary  system  displaced  in  the  negative  direction  0.1 


Initial  phase  angle*  .0000000000+00 

Initial  maximum  deflection*  - . 100000001 e+00 

natural  frequency*  . 250000000e+02  the  starting  time*  . 000000000e+00 
moment  of  Inertia  of  the  primary  mass*  .  1  09300001 e+00 

magnitude  of  the  secondary  mass*  .2 00000009e-02  length  of  the  primary  system* 
2.23000002  the  coefficient  of  restitution*  . 500000000e+00 

secondary  masses  Initial  velocity* 

.0000000000+00  the  gap  setting*  .2500000003+00 
desired  number  of  Impacts*  50 


mass*  . 200000009e-02 
effective  d=  .25 0000000e+00 
e*  . 500000000e+00 
x ( 0 ) *  -. 90000004 9e-01 
deltat*  . 999999978e-02 


the  time  Iteration  did  not  converge  for  I*  45 
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.0776 
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Impact*  5  errmom*  -.1117587090-07  errvel*  - . 2334  1 8579e-06 
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Impact*  10  errmom*  .1064739100-07  errvel*  -.1192092900-06 
Impact*  11  errmom*  - . 29C023224e-07  errvel*  . 2334 1 8579e-06 
Impact*  12  errmom*  - . 1 1 1 75C7u9e-07  errvel*  . 2364 1 8579o-0G 
Impact*  13  errmom*  -.5587930450-08  errvel*  . 2384 1 C579e-06 
Impact*  14  errmom*  - . 1 1 1 75G709e-07  errvel*  - . 2384 1 8579e-06 
Impact*  15  errmom*  .000000000e+00  errvel*  . 476837 1 58e-06 
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Fig.  4.  Output  from  Program  IDEAL 
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Fig.  4 — Continued 
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radians,  and  with  the  impacting  mass  in  contact  with  the 
positive  stop.  In  order  to  determine  how  important  these 
initial  conditions  were,  both  the  initial  position  and 
velocity  of  the  impacting  mass  was  varied.  Varying  the  ini¬ 
tial  position  of  the  impacting  mass  through  a  range  of  +1" 
to  -1"  from  the  center  of  the  gaps,  and  giving  it  a  velocity 
ranging  from  +  2  ft/sec  to  -2  ft/sec  led  to  a  ±4%  variation 
of  s  from  its  average  value.  Doubling  the  initial  ampli¬ 
tude  of  the  system  from  9  =0.08  radians  to  0  =0.16 

mQ  mQ 

radians  increased  s  by  20  percent.  None  of  these  changes 


significantly  affected  0  .  These  results  were  for 

m 

2  r 

I  =  0.1893  slug  -  ft  ,  impact  mass  =  0.001  slug,  L  =  2.23', 
d  =  3",  and  w  =  25  rad/ sec. 

initial  conditions  obviously  have  an  effect  on 
the  damper  efficiency  s.  This  effect  does  not  justify 
complicating  the  analysis  of  s  by  considering  initial  condi¬ 
tions,  especially  if  s  is  evaluated  while  keeping 
constant.  However,  the  variation  of  s  with  initial  ampli¬ 
tude  0^  suggests  that  the  system's  decline  in  amplitude 
is  not  perfectly  linear,  only  approximately  so. 

An  important  energy  consideration  for  space  opera¬ 
tions  is  that  the  only  energy  dissipation  will  be  due  to 
the  impacts.  Whatever  kinetic  or  potential  energy  is  lost 
due  to  these  impacts  will  be  converted  to  heat,  which  will 
be  distributed  between  the  primary  system  and  the  impacting 
mass.  This  heat  can  only  be  dissipated  through  radiation. 


which  is  a  very  slow  process.  Consider  a  primary  system 
with  inertia  I  about  its  rotation  point  and  a  natural  fre¬ 
quency  w  .  If  this  system's  only  energy  loss  is  due  to 
impact  damping,  which  reduces  the  system's  amplitude  from 
6m  to  0m  ,  then  the  change  in  total  energy  equals  the 
change  in  maximum  kinetic  energy,  so: 

iE  -  4T«x  -  iI(0O  )2  •  iI(®£  »2  (14) 

max  max 


where  0_  is  the  initial  maximum  angular  velocity,  and 
.  max 

0_  is  the  final  maximum  angular  velocity.  From  this  it 
°max 

is  easy  to  obtain: 

AE  =  ilu)  2(  6  2  -  6  2)  (15) 

2  n  '  mQ  mr  1  ' 


This  Ae  is  the  energy  converted  to  heat. 

A  worst  case  example  is  worked  out  for  the  labora¬ 
tory  model's  values  of  m,  I  and  w  ,  assuming  all  of  the 
heat  goes  to  raise  the  temperature  of  the  impacting  mass. 

For  the  smallest  mass  used,  m  =  0.000481  slug  (about  1/4 

2 

ounce),  I  =  0.188  slug  •  ft  ,  w  =  25.9  radians  per  second, 

and  9  and  0  are  chosen  to  be  0.10  radians  and  zero, 
mo  mr 

respectively.  This  gives: 


/  ^ 


AE  =  |(0.188) (25.9) 2(0.10) 2 


=  0.631  lb  •  ft 


=  8.103  x  10  BTU 


The  temperature  of  the  impacting  mass  will  then  be  raised 


AT°  = 


where  c  is  the  specific  heat  of  the  material.  If  c  = 
p  p 

0.109  BTU/ (lb  •  °F)  (the  value  for  nickel  steel  at  room 
temperature  according  to  Reference  24) ,  then 


AT°  =  0.048  °F 


While  this  temperature  increase  is  of  no  significance,  the 
fact  that  the  impacts  convert  kinetic  energy  to  heat  should 
be  remembered  when  designing  an  impact  damper,  especially 
if  a  very  small  mass  is  expected  to  absorb  a  great  deal  of 
energy . 


Forced  Vibration 

The  forced  vibration  impact  damper  depicted  in 
Figure  3  consists  of  an  upright  slender  beam  with  a  damper 
assembly  on  top  and  an  impacting  mass  free  to  move  between 
the  stops  of  the  damper  assembly.  A  time  dependent  angular 
displacement  is  applied  as  a  boundary  condition  to  the 
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bottom  of  the  beam.  The  motion  of  the  primary  system,  con¬ 
sisting  of  the  beam  and  damper  assembly,  was  not  obtained. 
However,  the  motion  of  a  similar  system  in  free  vibration, 
shown  in  Figure  5,  was  obtained  in  Reference  (25).  From 
this  solution,  the  natural  frequencies  can  be  found. 

The  undamped,  free  vibration  of  the  system  depicted 
in  Figure  5  is,  according  to  Reference  (25) : 


x(y,  t)  =  X (y)  e 


iu)_t 

F 


where : 


and: 


X (y)  =  C1  sin  (By)  +  C2  cos  (By) 


+  C3  sinh  (By)  +  C4  cosh  (By) 


“F  P 


sin  BL,  +  sinh  BL, 

c  =  _c  =  c  - 2 - k. 

2  4  ^lcos  BL,  +  cosh  BL, 

b  b 


C3  "  "C1 


is  determined  by  the  initial  conditions,  while  B  is 
solved  for  using: 


(cos  BL.  )  (cosh  BL,  )  +  Bp  (tanh  eLb 
b  b 


tan  BL,  )  =  0 
b 
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Mass,  modeled  as 
point  mass  at 
y=L 


Beam,  with 
mass/unit 
length  =  p 


y  coordinate,  moves 
with  _c enter  line 
of  beam 


x  coordin- 


7/7/7777777 


^  ate. 


fixed  in 
space 


Fig.  5.  Model  of  Free  Vibration  Problem 
Solved  in  Reference  25 


p  is  the  mass  per  unit  length  along  the  flexible  beam,  M 
is  the  total  mass  of  the  damper  assembly  at  the  top  of  the 
beam,  and  to  is  the  natural  frequency  of  this  system.  Once 
3  is  obtained  from  Equation  (20),  can  be  found  using 
Equation  (19)  . 


L'  . 


25 


III.  Laboratory  Models 
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Two  impact  damping  models  were  designed  and  con¬ 
structed  to  experimentally  study  the  performance  of  the 
damper  in  forced  and  free  vibration.  The  equations  of 
Chapter  II  were  derived  to  apply  to  the  models  depicted  in 
Figures  6  and  7.  The  free  vibration  model  of  Figure  6  was 
used  to  verify  the  analysis  of  Chapter  II.  The  forced 
vibration  model  of  Figure  7  was  used  to  take  measurements 
and  make  observations  on  its  motion.  Details  on  both  of 
these  models,  the  measurement  equipment  and  techniques,  and 
the  conversion  of  the  measurements  to  actual  displacements 
are  given  in  Appendix  B. 

Free  Vibration  Model 

The  laboratory  model  of  Figure  6  consisted  of  an 
aluminum  beam  suspended  by  a  short,  flexible  piece  of  steel 
acting  as  a  flex  plate,  with  the  damper  assembly  mounted  on 
the  bottom.  The  impacting  mass  is  suspended  from  a  point 
above  the  entire  system  to  minimize  friction  forces.  The 
quantities  needed  to  evaluate  the  motion  using  Equations 
(3)  and  (7)  are: 

I  =  0.188  slug  •  ft2 
c  =  0.02  lb  •  ft  •  sec 
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Flex  plate 


Aluminum  Beam 


The  model  used  to 
study  the  impact 
damper  in  free  vib) 
tion  (not  to  scale] 
The  top  of  the  fie) 
plate  was  attached 
a  secure  mount/  wh: 
the  rest  of  assemb! 
hung  underneath. 


Fig.  6.  Laboratory  Model  of  Free  Vibratio 
Impact  Damper 
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Quantities  needed  to  obtain  the  initial  conditions  after 
each  impact  using  Equations  (10)  through  (13)  are: 

L  =  2.21' 

e  =  0.4  -  0.5 

The  measurement  of  these  quantities  is  discussed 
in  Appendix  B,  but  some  comments  are  in  order  here.  The 
values  of  I,  WT,  r^,  dg,  I's,  m»  Lm  an^  L  were  measured, 
weighed,  and  calculated  to  an  acceptable  degree  of  accuracy 
with  little  uncertainty.  E  was  obtained  from  Reference  26. 
c  and  cm  were  obtained  by  measuring  the  reduction  in  ampli¬ 
tude  of  the  freely  oscillating  primary  system  and  impact¬ 
ing  mass  after  a  known  number  of  cycles;  c  and  c  could  then 

m 

be  calculated  by  the  logarithmic  decrement  method.  The 
value  of  c  obtained  in  this  manner  varied  from  0.01  to  0.04 
lbs  •  sec/ft,  c  =  0.02  was  taken  as  the  approximate  value. 

The  value  obtained  for  cm  did  not  vary  significantly  for 
different  tests,  but  the  impacting  mass  was  traveling  much 
slower  when  these  measurements  were  made  than  when  it  was 
used  in  the  impact  damper.  Since  damping  forces  are  not 
always  directly  proportional  to  velocity,  as  Equation  (6) 
treats  them,  this  could  be  a  source  of  error,  e  was  obtained 
by  allowing  each  of  the  impacting  masses  to  swing  as  a  pen¬ 
dulum  a  known  distance  and  strike  the  impacting  surface  of 
the  primary  system,  and  then  measuring  the  recoil  of  the 
primary  system  and  the  impacting  mass.  These  known 


quantities  and  measurements  were  converted  into  velocities 
for  the  primary  system  and  impacting  mass  both  before  and 
after  the  impact,  from  which  e  was  determined.  0.4  and  0.5 
are  the  upper  and  lower  values  of  e  obtained.  Finally, 
while  q  =  26.34  rad/ sec  was  the  calculated  value  for  the 
damped  frequency  of  the  system,  the  actual  frequency  of 
the  system  was  measured  as  25.9  rad/sec. 

The  motion  of  the  system  was  measured  using  two 
SR-4  type  AD-7  strain  gages,  centered  on  opposite  sides 
of  the  flex  plate.  These  gages  were  connected  to  a  Q-amp 
in  a  Type  535A  oscilloscope.  The  strain  e  obtained  in 
this  manner  could  be  converted  to  radians  of  displacement 
of  the  primary  system  using: 

6  =  4.80  x  10"5e  (21) 

where  e  is  measured  in  micro  inches  per  inch.  This  rela¬ 
tion  is  obtained  in  Appendix  B,  and  is  only  valid  for  small 
values  of  0,  where  the  relation  between  e  and  0  is  linear. 

Photographs  of  the  oscilloscope  trace  were  made  to 
obtain  the  beam  deflection  as  a  function  of  time.  Some 
of  these  photographs  are  shown  in  Figure  8 .  The  actual 
distance  between  the  impacting  surfaces  was  3"  in  this 
figure,  but  the  actual  distance  the  impacting  mass  could 
travel  between  the  impacting  surfaces  was  3"  minus  the 
diameter  of  the  impacting  mass.  The  diameter  of  the 
0.000481  slug  mass  was  15/32”,  the  0.00149  slug  mass  was 
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>; 
»  " 

a. 

* 


11/16",  and  the  0.00503  slug  mass  was  1-1/32".  d  is  the 
actual  horizontal  distance  between  impacting  surfaces  that 
the  impacting  mass  can  travel. 

It  can  be  seen  from  Figure  8  that  when  the  impact 
damper  is  operating,  the  maximum  angular  displacement 
reached  by  the  primary  system  decreases  linearly  with  time. 
This  is  not  the  constant  damping  of  an  exponentially  decay¬ 
ing  system.  The  usual  measures  of  damping,  such  as  log¬ 
arithmic  decrement,  will  not  be  independent  of  the  ampli¬ 
tude  of  the  system.  It  is  also  apparent  that  the  impact 
damper  ceases  to  be  effective  after  the  maximum  amplitude 
obtained  during  a  cycle  decreases  to  a  certain  value. 

This  value  is  denoted  0  ,  and  is  referred  to  as  the 

mr 

residual  amplitude.  When  the  damper  is  effective,  0^  > 

0  ,  where  0  is  the  maximum  amplitude  attained  by  the 

m^.  nij 

primary  system  during  cycle  j.  This  implies  that  the 
maximum  angular  displacement  attained  during  the  jth  cycle 
can  be  given  by: 


0  =  r  -  st . 

mj  3 


(22) 


where  9^  is  the  maximum  angular  deflection  attained  during 
cycle  j;  tj  is  the  time  at  which  the  system  reaches  maxi¬ 
mum  amplitude  during  cycle  j;  r  would  be  the  amplitude  0 


ITlr 


of  the  system  if  it  were  started  with  zero  velocity  at 


t.  =  tn  =0;  and  s  is  the  rate  at  which  the 
3  0 
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system's  maximum  amplitude  decreases  with  time,  which  will 
be  referred  to  as  the  damper  efficiency. 

Forced  Vibration  Model 

The  laboratory  model  of  a  continuous  system  in 
forced  vibration,  depicted  in  Figure  7,  consisted  of  a 
flexible,  upright  steel  beam,  with  a  time  dependent  angular 
displacement  applied  to  the  base,  and  a  damper  assembly  on 
top.  The  quantities  needed  to  obtain  the  natural  frequency 
from  equations  (16)  and  (17)  are: 


L,  =  18.5" 
b 

M  =  0.01382  slug 
p  =  0.00110  slug/in 
E  =  28  x  106  lb  •  in2 


1  .  4 

6144  ln 


For  these  values,  Equation  (17)  gives  thu  first  two  values 
of  6L  being  1.345  and  4.071.  Using  this,  Equation  (16) 
gives  the  first  two  natural  frequencies  as  being  37.28 
and  341.7  rad/ sec. 

A  model  MB  303  accelerometer  was  mounted  1"  below 
the  damper  assembly.  The  accelerometer  signal  was  ampli¬ 
fied  using  a  model  2614b  amplifier  powered  by  an  Endevco 
Mode  2621  power  supply;  this  signal  was  then  recorded  using 


a  Honeywell  Model  2106  visicorder.  High  frequency  noise 
required  sending  the  amplified  signal  through  a  simple 
low  pass  filter  before  it  reached  the  visicorder.  Details 
of  this  filter  are  given  in  Appendix  C. 
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IV.  Correlation  of  Analytical  and 
Laboratory  Results 

This  chapter  compares  the  computed  motion  of  the 
free  vibration  impact  damper  to  the  motion  measured  in  the 
laboratory.  Possible  sources  of  errors  in  both  laboratory 
measurements  and  in  the  attempt  to  preduct  the  motion  of  a 
real  system  using  LABSIM  and  the  equations  of  Chapter  III 
are  discussed. 

The  program  LABSIM,  explained  in  Chapter  II,  was 
used  with  the  measured  physical  quantities  of  the  labora¬ 
tory  free  vibration  model,  given  in  Chapter  III.  Table  1 

compares  values  of  s  and  0  obtained  from  LABSIM  for 

mr 

e  =  0.4  and  0.5  with  the  measured  laboratory  values  of  s 
and  0m  for  the  different  gap  settings  and  impacting  masses 
0.4  and  0.5  were  the  minimum  and  maximum  values  measured 
for  the  coefficient  of  restitution  e.  As  can  be  seen  the 
measured  value  of  damper  efficiency  s  is  never  more  than 
14  percent  greater  than  the  largest  value  of  s  computed, 
or  5  percent  less  than  the  smallest  value  computed.  How¬ 
ever,  the  measured  s  is  not  consistent  in  falling  between, 
above,  or  below  the  computed  values  of  s.  The  measured 
value  of  0^  may  be  more  than  twice  the  nearest  computed 
value.  The  rest  of  this  chapter  considers  possible  reasons 


for  these  imperfect  correlations. 


TABLE  1 


COMPARISON  OF  COMPUTER  SIMULATION  AND  LABORATORY 
RESULTS  FOR  s  and  6 


... 

Computer 

Laboratory 

mass 

(slugs) 

' 

B 

m 

e  = 

0.5 

d 

s 

6m 

_ r _ 

s 

6m 

r 

s 

CD 

0.000481 

3.53" 

0.0115 

0.0102 

0.0107 

0.0113 

0.010 

0.00149 

1.312" 

0.0134 

0.0040 

0.0110 

0.0032 

0.0152 

0.0082 

2.312" 

0.0199 

0.0071 

0.0176 

0.0048 

0.021 

0.0137 

3.312" 

0.027 

0.0105 

0.023 

0.0075 

0.028 

0.0184 

0.00503 

0.969" 

0.029 

0.025 

0.0016 

0.024 

0.0042 

1.969" 

0.043 

0.0033 

0.037 

0.0033 

0.038 

0.0063 

2.969" 

0.061 

0.0057 

0.052 

0.0026 

0.056 

0.0130 

Notes 

s  measured  in  radians/ sec. 

0  measured  in  radians. 
mr 

The  actual  distances  between  impacting  surfaces  in 
the  laboratory  were  set  at  2",  3",  and  4";  the  values  given 
here  for  the  gap  d  are  these  distances  minus  the  impacting 
mass's  diameter . 


Sources  of  errors  in  obtaining  the  laboratory 

results  fall  into  two  categories:  (1)  errors  in  measuring 

the  physical  parameters  of  the  system,  and  (2)  errors  in 

measuring  the  resulting  output.  Of  the  physical  parameters 

measured,  c  and  c  are  the  most  uncertain,  for  the  reasons 

m 

discussed  in  Chapter  III.  In  measuring  the  mass  of  the 
impacting  masses,  the  0.00005  slug  mass  of  the  line  support 
ing  them  was  neglected.  This  certainly  increased  the  effec 
tive  mass  of  the  impacting  masses  by  some  small  amount. 
Measured  lengths  under  six  inches  could  have  up  to  1/32" 
error  in  them,  measured  lengths  over  six  feet  could  have 
up  to  one  inch  error  in  them.  In  calculating  I,  the  dis¬ 
tributed  mass  of  the  damper  assembly  was  modeled  as  a 
point  mass  26"  from  the  rotation  point.  These  are  only 
some  of  tha  error  sources,  most  of  which  can  be  assumed 
to  be  small.  With  the  exception  of  c,  c^,  and  e,  a  quali¬ 
tative  estimate  of  the  errors  in  the  values  of  the  physical 
parameters  given  in  Chapter  III  would  be  that  errors  are 
±5  percent  of  the  quantities  given,  or  less. 

Some  of  the  error  sources  involved  in  measuring 
the  results  of  a  laboratory  test  are  simp’e:  the  strain 
gage  used  was  accurate  to  within  ±2  percent,  while  the 
traveling  microscope  used  in  measuring  the  photographed 
oscilloscope  trace  had  a  small  amount  of  play  in  the 
adjustment,  causing  errors  of  approximately  0.1  percent 
or  less.  Human  judgement  provided  another  error  source; 
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in  particular,  in  measuring  the  photographed  oscilloscope 

trace  one  had  to  decide  where  to  take  a  measurement  from  on 

an  often  fuzzy  trace  edge.  A  qualitative  estimate  of  the 

errors  in  measuring  s  and  0  is  that  the  smallest  values  of 

r 

s  and  0  may  be  in  error  by  as  much  as  ±20%,  with  most 
r 

values  of  s  and  0  being  accurate  to  within  ±5%  or  less. 

m 

r 

The  computer  model  of  the  impact  damper,  LABSIM, 

gives  the  correct  values  for  s  and  0  for  the  numbers  it 

mr 

is  given  and  the  operations  it  performs.  Errors  caused 

by  limitations  in  the  accuracy  of  single  precision  FORTRAN 

would  be  insignificant  (less  than  1  percent  error)  for  the 

numbers  and  operations  employed.  The  only  reason  LABSIM 

would  not  give  the  motion  of  the  laboratory  model  of  the 

free  vibration  impact  damper  would  be  if  the  equations  of 

motion,  or  their  solutions,  for  this  system  are  in  error. 

The  equations  of  motion  of  the  primary  system 

(Equation  2)  and  the  impacting  mass  (Equation  6)  employed 

in  LABSIM  are  based  upon  the  assumption  that  restoring 

moments  and  forces  (k0  and  k^x)  are  proportional  to  dis- 

•  • 

placement,  and  damping  moments  and  forces  (c0  and  c^x) 
are  proportional  to  velocities.  It  is  commonly  accepted 
that  spring  forces  on  a  body  and  gravity  forces  on  a  pendu¬ 
lum  restricted  to  small  displacements  are  both  approxi¬ 
mately  proportional  to  displacement.  Damping  forces  are 
not  as  well  understood,  and  resistance  to  motion  is  often 
assumed  to  be  independent  of  all  but  the  direction  of 


motion  (dry  friction)  or  to  be  proportional  to  velocitv 
squared  (aerodynamic  drag) .  Thus,  the  assumed  damping 
forces  in  Equations  (2)  and  (6)  may  differ  from  the  actual 
damping  of  the  laboratory  model. 

A  certain  source  of  error  in  the  attempt  to  model 
a  real  system  using  LABSIM  lies  in  the  small  angle  approxi¬ 
mations  made  in  the  derivation  of  Equations  (2)  and  (6) . 
Other  differences  between  the  LABSIM  model  and  reality  are 
the  assumption  that  impacts  are  of  infinitely  short  dura¬ 
tion,  or  that  an  impact  occurs  whenever  an  iteration  puts 
the  impacting  mass  within  0.000001  feet  of  a  stop.  Finally, 
it  is  unlikely  that  Equations  (2)  and  (6)  take  into 
account  all  of  the  forces  acting  upon  the  primary  system 
V  and  the  impacting  mass.  It  would  be  very  difficult  to 

quantify  all  of  these  sources  of  error,  or  to  say  if  these 
errors  add  up  or  cancel  out  over  many  cycles. 

In  view  of  the  errors  mentioned,  the  agreement 
between  the  computed  and  measured  values  of  s  seems  accept¬ 
able.  Which  of  these  errors  causes  the  differences  between 

the  computed  and  measured  values  of  0  is  unknown.  While 

r 

LABSIM  does  a  poor  job  of  predicting  0  ,  it  is  reasonable 

r 

to  assume  that  the  errors  in  LABSIM  do  not  significantly 
favor  one  set  of  system  parameters  over  another.  There¬ 
fore,  for  an  impact  damper  with  equations  of  motion  similar 
to  Equations  (2)  and  (6),  LABSIM  should  be  able  to  predict 
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how  changes  in  system  parameters  will  affect  s  and  0  , 

mr 

even  if  it  does  not  reach  a  correct  actual  value  of 


V.  Results  and  Discussion 


Free  Vibration 

All  of  the  results  discussed  here  were  obtained 

by  using  the  computer  model  of  the  impact  damper  IDEAL, 

the  model  with  no  gravity  or  viscous  damping.  The  system 

modeled  had  the  parameters: 

0  =  initial  amplitude  =  0.10  radian; 

m0  2 
I  =  moment  of  inertia  =  0.1893  slug  •  ft  ; 

L  =  impact  moment  arm  =  2.23  ft; 

d  =  gap  setting  =  0.25  ft; 

m  =  impact  mass  =  0.001  slug; 

e  =  coefficient  of  restitution  =  0.5;  and 

ojn  =  natural  frequency  =  25  radians/sec. 

The  parameters  0  and  I  were  kept  constant,  as  was  the 

m0 

flex  -  plate  stiffness;  all  of  the  other  parameters  were 
varied  one  at  a  time.  The  effect  that  changes  in  the  last 
five  parameters  had  on  the  damper  efficiency,  s,  and  the 
residual  amplitude,  0^  are  shown  graphically. 

Figure  9,  which  plots  damper  efficiency  s  for  dif¬ 
ferent  impact  moment  arms  L,  was  made  by  varying  L  from 
1  ft  to  3-1/4  ft,  in  quarter  foot  increments,  while  holding 
all  other  parameters  at  the  values  given  in  the  preceding 
paragraph.  Figure  10,  which  gives  s  versus  d,  was  obtained 


rad/ sec 
016-1- 


similarly,  by  varying  the  gap  setting  d  from  1  inch  to 
4  inches  in  quarter  inch  increments,  while  all  other 
parameters  were  held  constant.  The  same  technique  was 
used  in  obtaining  Figure  11,  s  versus  coefficient  of  restitu 
tion  e,  and  Figure  12,  s  versus  the  natural  frequency  w  . 

The  data  for  Figure  13  was  obtained  in  the  same  manner  as 
for  the  parameters  of  Figures  9  to  12,  but  the  plot  was 
made  a  little  different.  It  was  found  that  the  damper  effi¬ 
ciency  s  is  approximately  proportional  to  the  quantity: 


I  +  mL 


The  data  for  Figure  13  was  obtained  by  varying  m  from  0.001 

slug  to  0.005  slug  in  increments  of  0.005  slug.  The  impact 

moment  arm  L  was  kept  at  2.23  ft  and  the  moment  of  inertia 

2 

I  was  kept  at  0.1893  slug  •  ft  . 

The  approximate  linearity  of  Figures  9,  10,  12,  and 

13  and  the  linearity  of  Figure  11  for  e  >_  0.3  make  the 

analysis  of  s  a  simple  matter.  Since  s  appears  to  be 

directly  proportional  to  the  impact  moment  arm  L,  the  gap, 

d,  one  minus  the  coefficient  of  restitution,  1  -  e  (for 

e  greater  than  or  equal  to  0.3),  the  natural  frequency  w  , 

2 

and  the  mass  over  the  total  inertia,  m/(I+mL  ),  the  follow¬ 
ing  relation  can  be  written  for  e  greater  than  or  equal  to 


fmLd(l-e) 1 


.  I+mL2 
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Fig.  12.  Damper  Efficiency  for  Different  Values  of  Natur. 
with  m  =  0.001  slug,  I  =  0.1893  slug  •  ft^,  L  =  2.23',  d  =  3",  an* 


It  is  interesting  to  point  out  that  the  damper 
efficiency  improves  as  the  coefficient  of  restitution  e 
decreases,  even  for  e  less  than  0.3,  as  Figure  11  shows. 
This  is  not  surprising,  since  low  coefficients  of  restitu¬ 
tion  result  in  a  greater  energy  loss.  However,  according 
to  References  (7),  (8),  (10),  (16),  and  (20),  for  an  impact 
damper  used  in  forced  vibration,  optimal  damper  efficiency 
results  by  choosing  e  as  close  to  unity  as  possible.  The 
difference  between  the  two  systems  is  that  for  the  forced 
vibration  in  steady  state  motion,  the  forcing  function 
must  provide  a  constant  energy  input  equal  to  the  energy 
loss,  through  impact  or  any  other  mechanism.  In  free 
vibration,  the  system  starts  with  a  certain  total  energy 
and  any  energy  lost  is  not  restored.  This  reduces  the  pos¬ 
sible  motion.  This  illustrates  the  fact  that  parameters 
given  in  the  literature  which  optimize  the  impact  damper's 
performance  in  forced  vibration  do  not,  in  general,  opti¬ 
mize  the  impact  damper's  performance  in  free  vibration. 

A  useful  feature  of  Equation  (23)  is  that  the  term 

in  brackets  is  dimensionless.  This  requires  that  cg  be 

dimensionless.  If  both  s  and  w  are  measured  in  radians 

n 

per  second,  the  solution  of  Equation  (23)  for  cg,  taken  as 
the  average  solution  for  cg  over  the  full  range  of  param¬ 
eter  variations,  is: 


(24) 
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n  v,  _  s  ( I+mL  ) 

U*'3J  mLd(l-e)u) 

n 

Stated  another  way,  the  damper  efficiency  s  is  related  to 
its  parameters  by  the  relation: 


s 


mLd(l-e)(jo 
_ n 

I+mL2 


(25) 


This  relation  is  only  good  for  e  21  0.3,  and  may  not  apply 
for  parameters  outside  the  range  of  those  used  here  to 
obtain  Equation  (25).  Within  these  constraints.  Equation 
(25)  should  allow  one  to  choose  a  damper  efficiency  s. 

Figures  14  through  18  are  plots  of  0  versus  L,  d, 

2  ^ 

e,  u)n,  and  m/(I+mL  ),  obtained  in  the  same  manner  as  Figures 
5  through  9.  These  graphs  illustrate  that  the  relationship 
of  the  residual  amplitude  0^  to  the  parameters  varied  is 
complicated.  Attempts  were  made  to  solve  for  0^  as  a  func¬ 
tion  of  L/d,  but  these  results  were  considerably  more 
erratic  than  those  of  Figure  14,  where  L  is  varied  and  d 
kept  constant;  and  Figure  15,  where  d  is  varied  and  L  kept 
constant . 

By  plotting  the  results  given  in  Figure  14  on 
logarithmic  graph  paper  (also  called  log-log  graph  paper) , 
the  graph  of  Figure  19  was  obtained.  This  graph  implies 
that  0^  can  be  approximated  as: 


Fig.  14.  Residual  Amplitude  for  Different  Impact  Moment  Arms,  with 
.001  slug,  I  =  0.1893  slug  •  ft^,  d  =  3”,  e  =  0.5,  and  u  =25  rad/sec 


rad 


Fig.  17.  Residual  Amplitude  for  Different  Values  of  Natural  Fre¬ 
quency,  with  m  =  0.001  slug,  I  =  0.1893  slug  •  ft2,  L  =  2.23',  d  =  3",  an 
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Fig.  19.  Residual  amplitude  for  Different  Values 
of  the  Impact  Moment  Arm  on  a  Base  10  Logarithmic  Scale, 
m  =  0.001  slug,  I  =  0.1893  slug  •  ft2,  d  =  3",  e  =  0.5, 
and  a)  =25  rad/sec 
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where  L,  measured  in  feet,  is  divided  by  one  foot  to  non¬ 
dimens  ionalize  it.  If  we  assume  that  as  L  approaches 

infinity,  8  approaches  zero,  or  some  negligible  value, 
mr  . 

this  relation  becomes: 


0.022 


(26) 


It  is  emphasized  here  that  this  is  an  approximation  to  a 
function  that  is  not  well  understood,  and  this  approxima¬ 
tion  is  obviously  invalid  for  small  values  of  L. 

Figure  15  shows  that,  with  the  exception  of  one 
data  point,  9^  is  almost  exactly  linear  in  the  gap  d. 

This  one  bad  data  point  can  easily  be  explained.  The 
amplitude  where  the  impact  damper  ceases  to  function  was 
not  always  precisely  defined  in  the  computer  simulation. 

In  many  cases  there  was  a  sharp  transition  from  where  the 
amplitude  declined  linearly  to  where  the  damping  action 
ceased.  However,  in  some  cases  the  damper  transition  from 
effective  to  ineffective  operation  took  place  over  one  or 
two  cycles,  making  the  determination  of  9^  something  of  a 
judgement  call.  For  this  reason,  occasional  variations 
from  what  appears  to  be  an  otherwise  well-defined  trend  can 
be  expected.  In  the  case  of  9^  versus  d,  it  can  be 
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stated  with  confidence  that  0  is  linear  in  d. 

mr 

Figure  16,  which  is  a  graph  of  residual  amplitude 

0  versus  coefficient  of  restitution  e,  clearly  shows  that 
mr 

a  high  value  of  e  is  useful  in  minimizing  0  .  Unfortun- 

r 

ately,  a  low  value  of  e  is  desired  to  maximize  damper 
efficiency  s.  A  quadratic  function  was  fit  to  the  points 
of  Figure  16,  using  the  least  squares  method  to  minimize 
errors.  This  quadratic  function  is: 

0  (e)  =  0.02679  -  0.05132e  +  0.0248e2  (27) 

mr 

This  curve  comes  very  close  to  all  of  the  points  plotted, 

but  it  is  only  an  approximation  to  an  unknown  function 

relating  0  to  e . 
mr 

The  plot  of  6  versus  natural  frequency  to  ,  shown 
m^,  n 

in  Figure  17,  shows  0m  to  be  unaffected  by  By 

examining  the  equations  of  Appendix  A  upon  which  the  com¬ 
puter  simulation  is  based,  it  is  seen  that  an  increase  in 
wn  causes  a  proportional  increase  in  the  rate  at  which  the 
system  operates,  but  does  not  otherwise  affect  how  it 
operates.  Therefore,  it  is  logical  that  damper  efficiency 

s  would  be  proportional  to  w  ,  but  that  0  would  be 

n  m 

r 

unaffected  by  w  . 

J  n 

In  evaluating  Figure  18,  the  plot  of  0  versus 

2  mr 
m/(I+mL  ),  it  is  difficult  to  envision  what  occurs  as  m 

2  2 

becomes  very  large,  driving  m/(I+mL  )  to  the  value  of  1/L  . 


different  value  for  m  probably  warrants  a  recalculation 
of  C q.  Cq  is  also  only  correct  if  L  and  d  are  measured 
in  feet,  though  it  could  easily  be  recalculated  for  other 
units . 

Forced  Vibration 

The  motion  of  the  beam  in  forced  vibration  was  mea¬ 
sured  with  an  acclerometer  mounted  just  below  the  damper 
assembly.  The  output  of  the  accelerometer  had  to  have  high 
frequency  noise  filtered  out  before  the  visicorder  trace 
of  this  output  would  become  readable.  This  noise  was  seri¬ 
ous  enough  to  drown  out  the  sinusoidal  signal  expected  when 
the  beam  was  excited  without  the  impact  damper.  The  noise 
problem  was  even  more  serious  when  the  impact  damper  was  in 
place.  Presumably,  the  accelerometer  was  picking  up  the 
high  frequency  beam  vibrations  that  caused  acoustical  noise 

A  hand  tracing  of  the  filtered  accelerometer  output 
is  shown  in  Figure  20,  for  the  motion  without  the  impacting 
mass,  and  in  Figure  21,  for  the  motion  with  a  0.00149  slug 
impacting  mass  and  a  total  gap  of  two  inches.  The  actual 
visicorder  traces  are  given  in  Appendix  B.  Figure  21  sug¬ 
gests  that  there  is  no  simple  steady-state  motion  present. 
This  suspicion  was  confirmed  by  the  sound  of  the  irregu¬ 
larly  occurring  impacts.  Other  than  this,  the  visicorder 
traces  were  of  little  use  in  a  quantitative  analysis  of  the 
forced  system. 


Fig.  20.  Trace  of  Visicorder  Output  for  Forced  Vibration  Model  with 
no  Impacting  Mass 


Simply  observing  and  listening  to  the  forced  sys¬ 
tem  with  the  damper  operating  gave  some  valuable  informa¬ 
tion.  The  motion  of  the  system  was  not  enough  to  initiate 
and  sustain  impacts  unless  the  system  was  forced  near  its 
first  resonant  frequency,  or  unless  the  gap  setting  was 
very  small.  Higher  resonant  frequencies  caused  very  low 
amplitude  vibrations.  No  attempt  was  made  to  judge  the 
damper's  effectiveness  at  these  resonant  frequencies. 

At  frequencies  near  the  first  resonance  and  gaps 
greater  than  one  inch  the  0.00149  slug  mass  did  reduce  the 
amplitude  of  the  system  by  a  factor  of  at  least  two,  but  it 
did  not  approach  any  detectable  steady  state  operation. 

For  very  small  gap  settings  there  was  a  possible  steady 
state  reached,  but  no  detectable  reduction  in  amplitude. 
Using  the  0.00503  slug  mass  at  any  gap  large  enough  for  it 
to  affect  the  motion  of  the  system  led  to  a  very  erratic 
motion  of  the  system  with  no  evident  steady  state  operating 
state,  and  no  significant  sustained  reduction  in  amplitude. 
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I.  Conclusions  and  Recommendations 


The  impact  damper  shows  considerable  potential  in 
reducing  the  free  vibration  of  long,  lightly  damped  struc¬ 
tures.  It  is  especially  promising  for  long  structures 
since  the  effectiveness  of  the  damper  increases,  and  the 
amplitude  at  which  the  damper  becomes  ineffective  decreases 
as  the  damper  is  moved  farther  from  the  rotation  point. 

For  a  coefficient  of  restitution  of  0.3  or  greater,  Equa¬ 
tion  (24)  gives  a  good  estimate  of  the  rate  at  which  the 
oscillations  will  be  reduced,  while  Equation  (29)  gives  an 
estimate  of  the  amplitude  where  the  damper  becomes  ineffec¬ 
tive.  Even  in  the  ineffective  region,  where  impacts  occur 
sporadically,  each  impact  converts  some  of  the  kinetic 
energy  of  the  system  to  heat,  so  only  after  all  impacts 
cease  does  the  damper  become  totally  ineffective. 

A  problem  with  the  impact  damper  is  the  impacting 
mass's  need  for  room  to  travel  and  stops  to  impact  against. 
If  this  cannot  be  designed  into  or  added  onto  the  structure 
without  an  unacceptable  gain  in  weight  or  loss  of  struc¬ 
tural  strength,  the  impact  damper  should  not  be  used.  If 
structural  strengtu  is  a  problem,  using  a  damper  with  a 
very  low  coefficient  of  restitution  might  be  a  solution. 
While  Equation  (6)  does  not  hold  for  e  less  than  0.3, 


Figure  11  indicates  that  damper  efficiency  would  still  be 
higher  than  for  any  e  greater  than  0.3.  However,  Figure  16 
indicates  that  the  residual  amplitude  also  will  be  high. 

The  residual  amplitude  of  the  damper  can  be  elimi¬ 
nated  using  other  damping  techniques,  or,  to  generalize  an 
idea  proposed  in  Reference  (23),  by  putting  two  dampers  in 
parallel.  One  damper  would  be  designed  to  quickly  reduce 
large  amplitude  oscillations,  while  the  other  would  con¬ 
tinue  to  reduce  oscillations  to  a  smaller  value  at  a  slower 
rate  after  the  fast  damper  became  ineffective.  Another  pos 
sibility  is  a  system  which  reduces  the  gap  as  the  ampli¬ 
tude  of  the  primary  system  decreases.  Maximum  damper 
efficiency  could  be  obtained  keeping  the  gap  as  large  as 
possible  without  the  damper  becoming  ineffective;  i.e., 
keep  the  gap  just  small  enough  so  impacts  are  sustained. 

Another  potentially  useful  variation  on  the  basic 
design  of  the  impact  damper  would  be  to  let  the  impacting 
mass  travel  in  two  or  three  dimensions,  impacting  against 
a  ring  or  sphere  enclosing  it.  This  could  be  of  use  in 
damping  out  oscillations  about  more  than  one  axis.  This 
same  problem  could  be  attacked  by  orienting  one-dimensional 
impact  dampers  along  all  possible  rotation  planes,  but  this 
could  get  into  weight  problems. 

Other  variations  on  the  impact  damper  in  free  vibra 
tion  would  be  to  replace  the  impacting  mass  with  many 
masses  or  a  liquid.  Or,  the  hard  stops  could  be  replaced 


by  springs  or  dashpots.  Many  of  these  variations  have  been 
studied  for  the  impact  damper  in  forced  vibration,  but,  as 
has  been  noted,  what  maximizes  the  efficiency  of  the  forced 
damper  does  not  necessarily  maximize  the  efficiency  of  the 
free  damper.  The  basic  system  studied  here  shows  enough 
potential  to  warrant  further  study. 

With  the  exception  of  Reference  (5),  all  references 
on  the  subject  agreed  that  the  impact  damper  shows  great 
promise  in  eliminating  forced  vibrations.  (Reference  (5) 
studied  impact  dampers  solely  for  the  purpose  of  elimi¬ 
nating  vibrations  in  ship's  hulls.)  The  results  of  this 
study  shows  that  the  impact  damper  can  provide  some  damping 
of  structures  in  forced  vibration  near  resonance.  However, 
this  damping  was  sensitive  to  system  parameters,  and  no 
steady-state  motion  was  found. 

The  forced  motion  of  the  impact  damper  needs 
further  study,  not  only  to  resolve  differences  in  theories, 
but  also  in  the  laboratory.  Laboratory  models  should  be 
designed  not  only  with  the  objective  of  simulating  struc¬ 
tures  of  interest,  but  also  with  knowledge  of  the  limits 
of  measuring  equipment.  Many  instruments  are  poorly 
equipped  to  handle  vibrations  of  30  to  35  radians  per 
second. 

In  summary,  the  impact  damper  did  not  show  itself 
to  be  effective  or  predictable  in  reducing  or  eliminating 
force  vibrations.  However,  the  impact  damper  was  both 


effective  and  predictable  in  reducing  the  free  vibration  of 
a  structure  to  a  certain  value.  A  comparatively  small  mass 
can  greatly  reduce  the  amplitude  of  a  much  larger  vibrating 
structure,  in  many  cases  only  taking  a  few  cycles  to  do  so. 
The  impact  damper's  results  in  reducing  free  vibrations  not 
only  warrants  further  study,  but  also  warrants  careful, 
cautious  consideration  for  use  in  current,  applicable  vibra 
tion  problems. 
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Derivation  of  Equations 


Introduction 

The  equations  of  motion  and  their  solutions  for  the 
free  vibration  of  a  simple  one  degree  of  freedom  system  are 
well  known.  The  laboratory  model,  depicted  in  Figure  22, 
used  to  study  the  affect  of  the  impact  damper  on  free  vibra¬ 
tion  is  described  by,  for  the  period  between  impacts,  .he 
equations  of  motion  solved  for  later  in  this  appendix.  The 
model  depicted  in  Figure  23  is  a  continuous  system  equiva¬ 
lent  to  a  flexible  beam  with  a  lumped  mass  attached  to  one 
end.  The  free  vibration  of  this  type  of  structure  was 
solved  in  Reference  (25) ,  and  is  given  in  Chapter  II  of 
this  thesis.  These  solutions  can  be  evaluated  if  the 
position  and  velocity  of  the  system  are  known  for  some 
specified  time  t,  and  are  used  in  the  computer  models  of 
impact  dampers  in  Appendix  C. 

Position  and  Velocity 
Relations 

It  is  assumed  that  an  impact  can  be  modeled  as 
being  of  infinitely  short  duration.  During  this  kind  of 
impact,  for  the  system  of  Figure  24,  the  position  9(t) 
of  the  beam,  and  x(t)  of  the  impacting  mass  remain  unchanged, 


String  supporting  the  impacting  mass 


\ 

Fig.  24.  Coordinates^  and  Angles  Used  in 
Deriving  the  Equations  of  Motion  of  the  Free  Vibra¬ 
tion  Impact  Damper 


while  their  velocities  are  changed  instantaneously  from 
0  (t .  ^  ^  )  and  x  ( t .  ^  ^  to  0  ( t .  ^  )  and  x  ( t .  ^  )  .  While  an 
actual  impact  is  of  finite  duration,  this  time  is  so  small 
compared  to  the  periods  of  the  beam  and  suspended  mass 
that  this  assumption  is  justified.  This  causes  half  of 
the  initial  conditions  to  fall  out  immediately. 
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x  ( t . 

l 


0(ti 


(-) 

(-) 


) 

) 


The  remaining  two  initial  conditions  can  be  obtained  by 
the  conservation  of  angular  momentum  and  the  velocity  con¬ 
dition: 


e(x(ti(-))  -  0(ti(_))L)  =  0(ti(+))L  -  x(ti(+)) 

(30) 

Conservation  of  angular  momentum  requires  that  in  the 
absence  of  external  torques  the  total  angular  momentum  of 
a  system  about  a  fixed  point  0  remains  constant.  There 
are  external  torques  on  this  system,  so  the  magnitude  of 
the  total  angular  momentum  about  point  0  is  time  dependent, 
or : 


Hq  -  H0(t) 

However,  for  an  impact  of  infinitely  short  duration,  there 
is  only  a  momentum  exchange  within  the  total  system,  so: 


5 


E 


E 


3] 


It  is  convenient  to  take  the  fixed  point  0  to  be  the  point 
about  which  the  beam  rotates.  This  point  is  approximately 
fixed  for  small  angular  deflections  0.  This  will  be  shown 
to  be  true  when  the  motion  of  the  beam  is  solved  later  in 
this  appendix.  Strictly  speaking,  if  the  suspended  mass 
is  taken  to  be  a  point  mass,  then  the  magnitude  of  its 
angular  momentum  about  0  is: 


ft 


H, 


=  r .  x  mV 

flm  m 


(32) 


where  r_  is  the  vector  from  point  0  to  the  mass,  and  V 
Om  ^  m 


is  the  vector  velocity  of  the  mass.  Both  rQm  and  vary 


with  the  angular  deflection  a  of  the  mass,  but  if  a  is 
kept  small,  then  r^  and  Vm  are  approximately  perpendicular 


For  an  infinitely  short  impact,  rn  will  remain  constant 

Om 


in  magnitude,  so: 


=  mV  rn 
Om  m  Om 


(33) 


Two  further  useful  substitutions  are  obtained  by  noting 
that  for  small  a's,  the  vector  V^  is  essentially  aligned 
with  the  x  axis,  and  the  value  of  r^m  is  essentially  con¬ 
stant;  so,  writing: 


V  =  x(t) 
m 


6 


and 


rn  =  L 
Om 


the  angular  momentum  of  the  mass  becomes: 


H0m  =  ”Dt<t>L 


The  substitution  for  r^m  was  made  to  make  later  equations 
more  readable . 

The  magnitude  of  the  angular  momentum  about  point  0 
of  a  beam  with  moment  of  inertia  I  is  simply: 


H0b  =  IS (t) 


Conserving  angular  momentum  across  impact  i  results  in: 

10  ( t±  ^  — ^  )  +  mx(ti(~,)L  =  I0(ti(  +  ))  +  mx(tj.(  +  ,)L 


The  velocity  relations  can  be  rewritten  as: 


0(t.(  +  ))  =  fr{x(t.(  +  ))  +e[x(t.(“))  -  0(t.(_))L]} 

1  Li  1  X  X 


x(tj,(  +  ))  =  0(ti(  +  ))L  -  e[x(ti(  -  0(ti(_))L] 


Using  the  Equation  (37)  in  Equation  (36)  and  rearranging 


gives : 
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x(t.(  +  ))  =  re  ( t .  (  ))I(l+e) 

1  I+mtT  1 

+  x(t.(_) ) (mL  -  77)  ]  (39) 

Similarly,  using  Equation  (38)  in  Equation  (36)  gives: 

0(t.(+))  =  ~  2  [9(ti(_))(I  -  mL2e) 

1  I+ml/  1 

+  x(ti(-)  )mL(l+e)  ]  (40) 

These  relations,  along  with: 

0(t.(+))  =  9(ti(_))  (41) 

x(ti(+))  =  x(t.(_))  (42) 

give  the  initial  conditions  for  the  motion  between  impacts 
i  and  i+1  in  terms  of  the  final  conditions  between  i-1  and  i 

Motion  of  the  Vibrating  Beam 

The  equation  of  motion  for  a  rigid  beam  rotating 
about  a  fixed  point  0  can  be  obtained  from  the  relation: 

Z  Moments  =16  (43) 

The  beam  in  question  is  suspended  by  a  flex  plate  of 

length  d  ,  shown  in  Figure  25.  This  spring  can  be  said 
s 

to  apply  a  moment  of  magnitude  M  at  the  top  of  the  beam. 

This  moment  can  be  obtained  from  an  angle  0  from  the  rest 
(undeformed)  position  of  the  spring  steel  using  the  relation 
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’ig.  25.  Free  Body  Diagram  and  Coordinates  used  in  Deriving  Equa 
Motion  of  the  Free  Vibration  Impact  Damper 


(44) 


d  x  d0_  _  M(y) 

^.,2  <v  ei' 


where  M(y)  is  the  moment  of  distance  y  along  the  spring, 
and  the  x  and  y  coordinates  have  been  reversed  from  Refer¬ 
ence  (26).  This  equation  assumes  small  angles  6,  and  so 
assumes  small  deflection  in  the  x  direction  of  the  beam. 

The  free  body  diagram  shown  in  Figure  25  shows  that 
the  beam  moment  is  approximately  independent  of  y,  so  inte¬ 
grating  Equation  (44)  once  gives: 


0  (y)  = 


EI  ' 
s 


and  integrating  once  again  gives 


x(y)  =  \  |f-r 
s 


At  y  =  d  ,  0 (y)  is  the  angular  displacement  of  the  primary 
s 

system.  Evaluating  0  and  x  at  y  =  dg  gives: 

Md 

e'dsl=  ST'  <47> 

s 


,  Md 

X(ds)=  2  EI  '  2  0ds 

s 


If  the  aluminum  beam  oscillates  about  a  fixed  point  0  a 


distance  rg  from  the  top  of  the  beam,  then  for  small  values 
of  8  : 
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r  0=  r  (sin  0)  =  x(d  )  =  0d  (49) 

S  S  S  2  S 

so  r  =  id  ,  and  the  beam  rotates  about  a  point  located  at 

S  ^ 

the  center  cf  the  spring's  undeformed  position. 

The  moment  M  applied  to  the  beam  by  the  spring  is 
s 

opposite  in  sign  to  the  moment  at  the  end  of  the  spring. 
Using  Equation  (47),  this  gives: 


M 


s 


El 

s 


•e<as> 

d 


s 


(50) 


0  =  0 (d  )  is  the  deflection  from  vertical  of  the  beam, 
s 

Gravity  also  causes  a  moment  : 


M  =  -W_r  sin  0  =  -  Wm  r  0 

g  T  c  .g .  T  c  .g 


(51) 


where  WT  is  the  tot  1.  weight  of  the  primary  system,  and 

r  is  the  distance  from  the  center  of  gravity  of  the 

c  • 

primary  system  to  the  rotation  point  0.  Also,  there  is  a 
moment  due  to  damping  that  resists  the  motion  of  the  sys¬ 
tem.  This  damping  moment  is  assumed  to  be  proportional 
to  the  angular  velocity,  so  it  can  be  defined  as: 


M,  =  c0 


(52) 


The  equation  of  motion  of  the  system  then  becomes: 


ie  +  c0  +  ke  =  o 


(54) 


or : 


where 


k  =  W_  r 

T  c  .g , 


El  ' 
s 


(55) 


By  Reference  (27),  the  solution  to  0 ( t )  can  be  written: 


( -  —  t) 

8(t)  =  e'  21  '  (A  sin  qt  +  B  cos  qt) 


where  A  and  B  are  constants  depending  on  initial  condi¬ 
tions,  and: 


q  = 


(56) 


Since  we  are  interested  in  the  motion  of  the  system  between 
impact  i  at  time  t^  and  impact  i+1  at  time  t^+^,  the  solu¬ 
tion  is  more  conveniently  written  in  terms  of  time  t=t..+At 
where  At  ranges  from  0  to  t^+^-t^.  In  this  case: 

0(At)  =  e^  21  (A  sin  qAt  +  B  cos  qAt)  (57) 

If  the  angular  position  and  velocity  of  the  primary  system 
and  the  velocity  of  the  impacting  mass  are  known  immedi¬ 
ately  before  impact  i  at  time  t^  ^  ,  then  0(t^  +  ^)  and 
0(t^  +  ^)  are  given  by  Equations  (40)  and  (41). 

Using  these,  the  constants  A  and  B  can  be  solved, 
and  the  position  and  velocity  of  the  system  at  time  At  as 
At  ranges  from  t^  to  t^+^  is: 
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6  (At)  =  e("  21  At) 


e*(t.(+))  e(t.(+))c 

q  +  21^ 


sin  qAt 


+  0 ( ^ )  cos  qAt 


(58) 


c  <“  fi  At) 
0  (At)  =  (-  ^)e 


6 (t . (+) )  0(t  (+))c 

- q -  +  - 21^ -  Sin  qAt 


+  0(t^  +  ^)  cos  qAt 


+  e 


(-  £  At) 


0(t  (  +  k  ,  6(ti<+))C 

0(ti  )  +  - 21 


COS  qAt 


-  q0 (ti^~^ )  sin  qAt 


(59) 


If  the  system  is  undamped  these  equations  become: 

^i<+’’  (+, 

0 ( At)  =  - - -  sin  qAt  +  0(t^1  ' )  cos  qAt  (60) 


0(At)  =  0(t^+^)  cos  qAt  -  q0(t^  +  ^)  sin  qAt  (61) 


Motion  of  the  Secondary  System 

The  motion  of  the  secondary  system,  which  in  the 
laboratory  was  simply  a  steel  sphere  suspended  by  nylon 
thread,  is  most  easily  obtained  by: 

E  Moments  =  I^a  (62) 

In  this  case,  the  inertia  i  is  simply: 

Im  =  mLm2  (63) 
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The  only  moments  are  from  the  a  component  of  gravity: 


M  =  -mgL  sin  a  =  -  mgL  a 
g  *  m  ^  m 


and  the  moment  due  to  damping: 


Md  =  - 


c  a 
a 


(64) 


(65) 


so: 


2"  * 

mL  a  =  -  c  a  -  mgL  a  (66) 

m  a  ’  m 

In  order  to  conveniently  go  from  treating  the  mass 
as  a  damped  pendulum  to  treating  the  mass  as  being  free 
from  all  forces  except  the  impacts,  it  is  helpful  to  note 
that: 


x  =  L_  sin  =  L  a  (6  7) 

mam 

•  •  • 

x  =  aL  cos  a=  aL  (68) 

m  m 

x  =  aL  cos  a  -  aL  sin  a  =  aL  (69) 

m  mm 


This  implies  that  for  small  values  of  a,  little  accuracy 
is  lost  by  assuming  all  motion  is  in  the  x  direction,  so 
Equation  (66)  can  be  written: 

•  • 

mL  x  =  -  c  — - mgx 

m  a  L_ 


or 

••  • 

mx  +  c  x  +  k  x  =  0  (70) 


where : 


'm  T  2 
Lm 


k  =  m 

m 

IQ 


As  was  done  for  the  primary  system,  x  can  be 
solved  for  according  to  Reference  (27),  giving: 


(-  At) 

x(t)  =  e  (A  sin  q  t  +  B  cos  at) 

m  m  m  m 


where 


%  ■  -  & 


and,  writing  this  in  terms  of  At,  which  ranges  from  time  t 
to  time  t^+1: 

(-  ~  At) 

x(At)  =  e  (A  sin  q  At  +  B  cos  q  At) 

m  m  m  m 


A  and  B  can  be  solved  here  using  Equations  (39)  and  (42) 
mm 

the  same  way  A  and  B  were  solved  earlier  for  the  beam. 
Doing  this,  the  position  and  velocity  of  the  impacting 
mass  becomes: 


(-  Itlf/Xlt,1*1)  X<t)<+,|c  , 

*<“>  -  *  q_  +  -255 -  Sln 

\  ^m  ^m  / 


+  x(t.  ^  )  cos  q 


;  At] 


c  <-  •==  At)  x(t  x  ( t .  )  c 

x(At)  =  -Je  2m  —a—  +  ~ 2mq~ -  sin  qmAt 

m  m 


+  x(t.  ' )  cos  q  At 
1  m 


_<2H  At)  T/  *  <  +  )  X(ti(  +  ))cml 

+  e  2 '  x(t.{  >)  +  - 2S -  cos  qmAt 


-  x(t.'  ')qm  sin  qmAt  (76) 


In  the  undamped  case  these  equations  become: 

x(ti(+))  {+) 

x  (At)  =  - — | -  sin  qmAt  +  x(ti'  ' )  cos  qmAt 


x ( At)  =  x(t.(+))  cos  q^At  -  x(ti(+))  qm  sin  qmAt 


and  if  the  motion  of  the  mass  depends  only  upon  the  impacts 
i.e.,  no  gravity  or  damping  forces: 

x ( At)  =  x(ti(+))  +  x(ti(+))At  (79) 


x(At)  =  x  ( t . 


)  =  constant 


(80) 


Appendix  B 
Laboratory  Models 


Introduction 

This  appendix  describes  in  detail  the  laboratory 
models  and  equipment  used.  The  conversions  from  measured 
quantities  to  actual  displacements  are  developed,  as  well 
as  the  methods  used  to  indirectly  measure  some  of  the  sys¬ 
tems  parameters. 

Free  Vibration  Model 

The  model  used  for  the  free  vibration  experiments 
is  depicted  in  Figure  22  with  the  important  dimensions, 
masses,  and  properties  given  in  Figure  26.  The  1/8"  x  1" 
x  4-1/2"  steel  beam  used  as  a  flex-plate  (note  that  only  3" 
of  its  length  was  free  to  bend  as  a  spring)  was  attached 
by  screws  to  a  support  depicted  in  Figure  27.  This  sup¬ 
port  was  bolted  onto  a  1-1/2"  x  40"  x  42"  steel  plate  which 
was  itself  bolted  to  12"  x  12"  I-beams  which  extended  up 
from  the  building's  foundation. 

6  2 

Using  Equation  (50),  with  E  =  28  x  10  lb  •  in 

4 

and  I  =  1/6144  in  ,  the  moment  applied  by  the  spring  onto 
the  aluminum  beam  was  calculated  as  15196  lb  •  in.  The 
moment  of  inertia  of  the  beam  and  damper  assembly  was  cal¬ 
culated  by  modeling  the  damper  assembly  as  a  point  mass 
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Fig.  26.  Free  Vibration  Impact  Damper  Used 
in  Laboratory  with  Important  Parameters 


26"  from  the  assumed  rotation  point;  this  gave  I  =  0.1882 
2 

slug  •  ft  .  The  moment  of  inertia  was  then  measured  by 
hanging  the  beam  and  damper  assembly,  minus  the  steel 
spring,  from  1.5"  of  nylon  fishing  line  and  timing  it 
through  a  number  of  cycles  as  it  swung  as  a  pendulum. 
Neglecting  any  damping,  the  equation  of  motion  of  this 
system  is: 

16  +  rM^gS  =  0  (81) 

where  r  is  the  distance  from  the  rotation  point  to  the 
center  of  gravity  of  the  assembly,  and  M^,  is  the  total 
mass  of  the  assembly.  From  this,  the  natural  frequency  of 
the  system  is: 


was  measured  as  4.45  rad/ sec  and  M  is  0.08983  slug,  r 

was  calculated  to  be  15.4545",  and  g  was  taken  as  32.174 
2 

ft/sec  .  I  can  then  be  solved  for  using: 

rr^g 

I  =  — 2  (83) 

oj 

n 

.  2 

This  gave  I  =  0.1880  slug  *  ft  .  This  value  of  I  was  used 
in  all  calculations. 

Using  w  =  A/I,  Equation  (55),  and  r  =  15.7648, 
n  c  •  cj  • 

the  natural  frequency  of  the  beam,  was  calculated  to  be 
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26.34  rad/sec.  When  using  a  damping  factor  c  =  0.02, 
which  will  be  justified  later  in  this  appendix,  the  damped 
frequency  was  essentially  the  same.  The  actual  frequency 
observed  was  measured  as  25.9  rad/ sec.  The  difference  was 
assumed  to  be  the  result  of  small  inaccuracies  in  the  mea¬ 
sured  quantities.  For  the  purpose  of  the  computer  simula¬ 
tion,  the  measured  values  of  inertia  and  frequency  were 
used. 

The  damper  assembly  consisted  of  an  aluminum  bar, 
1/4"  x  1"  x  6",  attached  to  the  bottom  of  the  1"  x  1"  x  24 
aluminum  beam.  The  stops  were  made  of  steel  and  could  be 
attached  anywhere  along  the  aluminum  bar,  and  were  mounted 
as  shown  in  Figure  28.  The  steel  balls  used  as  impacting 
masses  were  hung  by  nylon  fishing  lines  attached  to  points 
84"  above  the  damper  assembly,  with  one  attachment  46" 
to  the  right  of  the  damper  assembly  and  the  other  50"  to 
the  left.  The  mass  was  hung  as  a  pendulum  to  minimize 
forces  other  than  impact. 

The  motion  of  the  beam  was  measured  with  two  SR-4 
Type  AD-7,  Lot  #B-32,  strain  gages  attached  to  the  steel 
spring,  centered  on  either  side.  These  strain  gages  had  a 
gage  factor  of  1.96  ±  2  percent.  The  strain  gages  were 
connected  to  a  Q-amp,  serial  number  002578  which  was 
installed  in  a  Type  535A  oscilloscope.  The  oscilloscope 
trace  was  photographed  using  an  oscilloscope  camera  C-12. 
The  peak-to-peak  amplitudes  on  the  photograph  were  then 
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Figure  28.  Details  of  Damper  Assembly 
and  Impacting  Surfaces 


measured  in  inches  by  a  traveling  microscope.  These  mea¬ 
surements  were  divided  by  the  measured  division  size  on  the 
photograph  to  give  the  amplitudes  in  scope  divisions.  The 
number  of  scope  divisions  was  multiplied  by  the  Q-amp 
setting  to  give  the  total  strain  of  the  two  strain  gages 
had  they  had  a  gage  factor  of  two.  Since  the  gage  factor 
was  1.96,  and  only  the  strain  on  one  side  of  the  stQel 
spring  was  desired,  this  measured  strain  was  co  srted 
to  the  actual  strain  e  using: 


e 


1,  2.0 

“  2 '1 . 96 


) 


£ 


m 


(84) 


This  strain  can  be  converted  into  radians  of  displacement 
using : 


e  =  - 


Mu 

EI_ 


(85) 


from  Reference  (26) ,  where  u  =  1/16"  is  the  distance  from 
the  neutral  surface  of  the  spring,  and  M  is  the  moment 
calculated  as  15199  lb  •  in.  This  gives 


e  =  0.020830 

£  =  (20,830  u'7")e  (86) 


So,  for  an  e  given  by  Equation  (86)  ,  the  angular  displace¬ 


ment  9  of  the  beam  would  be: 


6  =  207830  =  (4*80  x  1°5)e  (87) 

where  e  is  given  in  micro-inches/inch,  and  0  is  in  radians 

Measurements  were  made  by  manually  deflecting  the 
beam,  or  primary  system,  until  the  oscilloscope  trace  was 
at  the  desired  position  on  the  screen.  The  oscilloscope 
was  then  set  to  make  only  one  sweep  when  triggered.  The 
sweep  rate  for  most  measurements  was  0.5  sec/div  or  0.2 
sec/div.  The  oscilloscope  camera  was  then  fastened  into 
position  and  the  lens  opened.  In  rapid  succession  the 
primary  system  was  released  and  the  oscilloscope  was 
triggered,  so  the  camera  photographed  one  oscilloscope 
trace.  The  amplitude  of  the  photographed  cycles  was  mea¬ 
sured  in  inches  using  a  traveling  microscope.  These 
amplitudes  were  converted  to  scope  divisions  by  dividing 
by  the  measured  division  width,  then  converted  to  strain 
by  multiplying  by  the  Q-amp  strain  setting.  Equation  (82) 
was  then  used  to  get  the  maximum  angular  displacement  per 
cycle,  s  was  obtained  from  these  displacements  using  a 
linear  least  squares  fit. 

The  natural  damping  of  the  beam  assembly  was  mea¬ 
sured  by  deflecting  the  beam  a  desired  amount  and  allowing 
it  to  vibrate  freely  with  no  impacting  mass  in  place.  The 
photographed  oscilloscope  trace  was  then  used  to  measure 
its  damping  using  the  log  decrement  method.  From  this, 
the  damping  factor  c  of  Equation  (54)  was  calculated  using 
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where 


(89) 


(90) 


While  the  value  of  c  obtained  this  way  was  always  small, 
it  was  not  constant.  On  a  given  day,  £  appeared  to  be  a 
linear  function  of  the  initial  displacement,  but  this 
linear  function  changed  from  day  to  day.  This  is  illus¬ 
trated  in  Figure  29.  For  this  reason,  £  =  0.002  was  taken 
as  giving  a  good  average  value  for  the  amplitude  at  which 
most  measurements  were  taken  with  the  impact  damper  opera¬ 
ting.  This  £  gave  c  =  0.02  lb  •  ft  •  sec.  This  value  of 
c  was  used  in  the  computer  simulation  of  the  laboratory 
model . 

In  order  to  measure  both  the  damping  on  the  impact¬ 
ing  mass,  and  the  coefficient  of  restitution  e  between  the 
mass  and  the  stops,  the  position  of  the  impacting  mass  had 
to  be  measured  without  interfering  with  its  motion.  This 
was  done  by  mounting  two  pieces  of  white  poster  board  8" 
to  the  right  of  the  beam  assembly,  facing  the  beam  assembly 
Horizontal  and  vertical  lines  were  drawn  at  1"  intervals 


across  the  poster  board.  Twenty-six  feet  to  the  left  of 
the  beam  assembly,  a  lamp  was  pointed  towards  the  assembly. 
With  all  other  room  lights  dimmed,  the  impacting  mass  cast 
a  sharp  shadow  upon  the  poster  board.  The  impacting  mass 
was  removed  from  between  the  stops,  and  one  of  the  stops 
was  placed  on  the  end  of  the  damper  beam,  facing  out.  By 
watching  the  impacting  mass's  shadow  on  the  poster  board, 
the  mass  could  be  released  from  a  known  position  and  strike 
the  stop  at  a  known  position.  Assuming  negligible  damping, 
the  velocity  of  the  mass  before  impact  can  be  calculated 
using : 

mgAh  =  -^mV2  (91) 

where  Ah  is  the  difference  between  the  mass's  height  at 
release  and  its  height  at  impact.  The  distance  that  the 
beam  travels  due  to  the  impact  can  be  determined  from  the 
oscilloscope  trace,  this  maximum  angular  deflection  times 
the  natural  frequency  of  the  beam  gave  its  maximum  angular 
velocity,  which  occurred  immediately  after  the  impact. 

The  velocity  of  the  impacting  mass  after  impact  was  calcu¬ 
lated  by  using  conservation  of  momentum,  Equation  (36) . 

The  coefficient  of  restitution,  then,  becomes: 


6(  +  )L  - 

e  =  -  V(_) 


(92) 
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where  (+)  implies  immediately  after  the  impact,  (-)  implies 


immediately  before  the  impact,  L  is  the  distance  from  the 

rotation  point  of  the  beam  to  the  impact  height,  and 
*  ( -) 

0  =0.  Using  this  method  the  coefficient  of  restitu¬ 

tion  for  the  steel  balls  striking  the  steel  stop  was  found 
to  be  between  0.40  and  0.50,  as  is  seen  in  Table  2. 


TABLE  2 

QUANTITIES  FOR  CALCULATION  OF  COEFFICIENT 
OF  RESTITUTION 


M 

Ah 

v<-> 

e 

0 

m 

0<+) 

v(+) 

e 

0.000481 

11-3/4 

7.94 

100 

0.00254 

0.0657 

-3.80 

0.50 

0.00149 

11-1/2 

7.85 

280 

0.00711 

0.1841 

-2.77 

0.40 

0.00503 

11-3/8 

7.81 

900 

0.0228 

0.592 

-2.31 

0.46 

Data  from  which  coefficient  of  restitution  e  is 

calculated,  m  is  in  slugs.  Ah  is  in  inches,  and 

are  in  feet/ sec,  e  is  in  micro-inches/ inch,  0  is  in  radians 
.  ,  +  v  m 

and  0  =  0  uj  is  in  radians/ sec. 

m  n 

The  same  lamp  and  poster  board  arrangement  was  used 
to  measure  the  damping  factor  on  the  impacting  mass,  but 
without  the  beam  assembly  in  place.  The  mass  was  dis¬ 
placed  to  a  known  position  and  released.  After  a  known 
number  of  oscillations,  its  maximum  amplitude  was  noted  and 
its  damping  c  was  solved  using: 


and  L  =  84"  is  the  difference  in  height  from  where  the  mass 
is  at  rest  and  the  points  from  which  it  is  suspended.  The 
damping  measured  in  this  way  had  an  amplitude  dependence; 

C  increased  with  amplitude.  The  values  of  c  were  calcu¬ 
lated  for  low  amplitudes  to  get  the  best  correlation  with 
the  position  of  the  mass  when  it  was  used  for  impact  damp¬ 
ing.  The  data  and  resulting  values  of  c  obtained  are 
shown  in  Table  3. 


TABLE  3 

QUANTITIES  FOR  CALCULATION  OF  VISCOUS 
DAMPING  FACTOR 


m 

YAi 

j 

y 

6 

c 

m 

0.00481 

4 

22 

0.0660 

0.0105 

0.0000217 

0.00149 

4 

61 

0.0231 

0.00368 

0.0000235 

0.00503 

4 

179 

0.00779 

0.00124 

0.0000267 

Data  from  which  the  viscous  damping  coefficient  of 
the  impacting  mass  is  calculated,  m  is  in  slugs,  and  c  is 
in  lb  •  sec/ft.  m 


Forced  Vibration  Model 

The  model  used  for  the  forced  vibration  experi¬ 
ments  is  depicted  in  Figure  30  with  the  important  dimen¬ 
sions,  masses,  and  properties  given  in  Figure  31.  The 
first  resonant  frequency  was  calculated  according  to  Refer¬ 
ence  (25) ,  where  the  motion  of  a  beam  in  free  vibration 
with  one  end  clamped  and  a  point  mass  at  the  other  end  was 
solved.  The  beam  was  mounted  upon  a  block  clamped  upon  the 
pin  of  a  pin-tree  beam  test  apparatus.  This  test  apparatus 
could  rotate  +24°  to  -24°  at  a  frequency  range  of  0  to 
83.3  Hz  (0  to  5000  RPM) .  For  computational  purposes  the 
beam  was  treated  as  if  its  end  extended  to  the  center  of 
the  pin,  where  a  known  sinusoidal  angular  displacement  was 
applied. 

The  damper  assembly  was  simply  an  aluminum  bar, 

1/4"  x  1"  x  5"  mounted  on  top  of  the  vibrating  beam,  with 
the  stops  mounted  as  shown  in  Figure  28 .  The  stops  were 
made  of  steel  and  could  be  attached  anywhere  on  the  alumi¬ 
num  bar.  The  steel  ball  used  as  the  impacting  mass  was  sus¬ 
pended  by  two  pieces  of  fishing  line  attached  to  points 
93"  above  and  74"  to  one  side  and  90"  above  and  75"  to  the 
other  side  of  the  damper  assembly.  The  impacting  mass  was 
suspended  in  order  to  minimize  friction  forces  on  it,  so 
its  motion  could  be  treated  as  resulting  entirely  from  its 
impacts  with  the  stops.  When  the  mass  was  between  the 
damper  stops  it  could  displace  from  its  rest  position  a 
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Pig.  31.  Laboratory  Model  of  Forced  Vibra 
tion  Impact  Damper  with  Important  Parameters 


few  inches  at  best,  so  its  velocity  due  to  its  pendulum 
motion  was  very  small  compared  to  the  velocity  imparted 
to  it  by  the  impacts . 

An  accelerometer,  Model  MB  303,  serial  #  149235, 
was  mounted  on  the  beam  1"  below  the  bottom  of  the  damper 
assembly.  The  acclerometer  signal  was  amplified  using  a 
model  2614B  amplifier  powered  by  an  Endevco  Model  2621 
power  supply.  High  frequency  noise  was  filtered  out 
using  the  low-pass  filter  of  Figure  32,  and  the  signal  was 
then  recorded  using  a  Honeywell  visicorder  oscillograph 
Model  2106  with  an  M-1000  galvanometer.  The  output  was 
also  used  with  a  universal  counter  timer.  Model  726C,  to 
accurately  determine  the  frequency  of  the  system.  The 
resulting  visicorder  output  for  the  forced  vibration 
model  without  and  with  the  impacting  mass  is  shown  in 
Figures  33  and  34,  respectively. 


Signal 

Voltage 


150  MF 


1 20 


Visicorder 


O 


-o 


Fig.  32.  Low-Pass  Filter  Used  to  Filter 
the  Amplified  Accelerometer  Output  Before  Inputting 
it  Into  the  Visicorder 
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vertical  axis  represents  approximately .8  ft/ sec2.;  gap  =  1.31 


The  following  pages  contain  two  FORTRAN  77  computer 
programs  which  solve  for  the  motion  of  both  the  primary 
system  and  the  impacting  mass  for  an  impact  damper  in  free 
vibration.  The  first  program  is  the  ideal  case,  in  which 
the  primary  system  is  undamped  except  for  the  impact  damp¬ 
ing,  and  the  motion  of  the  impacting  mass  is  due  entirely 
to  the  impacts.  The  second  program  assumes  a  lightly 
damped  impacting  mass  hanging  as  a  pendulum.  Though  both 
programs  were  written  with  the  laboratory  model  in  mind, 
they  are  applicable  to  any  one  degree  of  freedom  system 
in  free  vibration  using  a  one  degree  of  freedom  impact 
damper . 

While  comments  explaining  the  programs  are  inserted 
in  appropriate  places,  a  few  additional  words  are  in  order. 
The  position  of  the  primary  system  was  put  in  the  form  of: 


u  (At)  =  a  sin  (mAt)  +  b  cos  (wAt) 


and 

-  2^  At 

(^(At)  =  e  [a  sin  (toAt)  +  b  cos  (u)At)  ] 


for  the  first  and  second  programs,  respectively.  The 
numbers  associated  with  parameters  in  these  programs 


were  assigned  with  units  of  feet,  seconds,  and  slugs 
intended . 

The  only  sources  of  error  in  the  program  are  the 
computer  round-off  errors,  and  the  errors  in  defining  an 
impact  as  occurring  whenever  the  impacting  mass  and  a  stop 
were  less  than  10  ^  ft  from  each  other.  This  leads  to 
small  errors  in  position  and  velocity  of  the  system  in 
impact.  In  order  to  judge  the  seriousness  of  these  errors, 
as  well  as  to  insure  the  equations  and  solution  approach 
are  correct,  the  errors  in  velocities  and  momentum  were 
calculated  after  each  impact,  using  Equations  (31)  and  (36) 
The  only  other  serious  source  of  errors  lies  in  the  approxi 
mations  and  assumptions  made  in  the  derivation  of  equations 
of  Appendix  A.  While  the  programs  will  generate  output 
for  any  magnitude  of  6^,  this  output  will  be  reasonably 
correct  only  if  the  assumption  of  small  angles  is  not 
violated . 
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program  IDEAL 

THIS  Pi' OCT  AM  FOLLOWS  THE  SYSTEM  THROUGH  A  SERIES  OF  IMPACTS 
IN  THE  "IDCAL"  CASE,  I.E.  WITH  NO  VISCOUS  DAMPING  OR 
CRAVI I Y  AFFECTS. 

Integer  I  ,  m ,  min  ,  c  ,  n  ,  J  ,  I  inper  r 

par  ai,i>.  tor  ( i»  =  SiJj  } 

real  phi  ( JC  :  n  )  .thetas!  1  :n)  ,  thetav!0:n)  ,  tJot  z  (  I  in)  , tdotv ( 0 : n )  , 
r  t  i  no  <  0  s  n )  , vo I  3 (  1  : n  *  ,  vo  1  v (0 : n  )  ,  cor.o  .  c  two  ,uli)rn  )  .b(flin)  , 

:  a  >  .  bb  ,  d theta  ,  1  thet  a  .  1  nc  t ,  1 1  I  me ,  t ,  L  t ,  t  ts  ,  at ,  tat ,  q  ,  S  . 

.  t-rrorUHn)  ,inaxorr  ,  var  I  ,  I  nip  t  (  1  :  n  )  , 

:  i  nor  t  ,  ua  s  u  ,  1  ,  e  ,  1  1  .  dci ,  qone  .  qtwo  , 

:  d  o  1 1  a  t  ,  II  .  I r , x ( 0 : n ) .vnllflin) , d . t  h  c  t  a  m  <  D : n 1 ,resld, 

i  er  r moil. ,  er  r  ve  1  ,  bniut  3  ,  dt  { 0  :  n  )  ,  mtt ,  mt ts  ,  m  i  dipt (  1  :  n  ) 

The  par amnt or s  varied  are  defined  below.  Only  one  set  of 
parameter-:  are  varied  hero.  In  practice,  "do  loops"  were 
used  to  obtain  a  largo  combination  of  parameters. 

mass-*0.002 
d  - 0 .  ?.  3 
c=0. 0 
qone  •-  Z  5 . 0 
1-2.  .13 

Conditions  from  any  previous  run  set  to  zero  in  the  loop  to  l.lne  101 

do  1  a  1  J ”  1  .  n . 1 
t no  t  •  ( J  1=0.0 

I  hoi  i  v I J  > " U . 0 
I  dot  .  <  J  )  =0 .0 
tdvt  >•  t  i  )  •  H .  0 
t  I  MO  j  )  !,i .  0 

vo 1  .(  i  )  n  .0 

vc  I  V  !  J  )  a .  0 

x  (  J  )  «.7.0 
their  Ml  j  )  -*0.0 
cent  i  i .  u  o 

T.ie  redlining  parameters  are  defined. 

.j t  |  1,1(0),  thetavl  0) ,  tlnieiO),  Inert, 

:  •  /  l  •/  1  ,  t-'otv!  0  /  / 

:  1  . 0 . -h .  1  .h.0.0. 1803, 

:  w . U . 0.0/ 

print*,  ’  ' 

print*,  *  ’ 

I  •  r  1  n  t.  ‘  .  ’ 

p  r  I  n  i  ,  ‘  I  n  1 1 1  -  1  phase  ar.g  1  e=  ’  ,  ph  1  I  0  1  ,  *  i  n  1 1 1  a  1 
:  m.j  v  i  mu-,;  *1'  i-jt  l*<ri-'.  t hot  v  '  0  >  .  '  nut ur ;  1  1'  r  oqurncy"  ’  ,  gone, 

I  '  t  h  •  .  Uil  l  '  li'J  t  Mil-  ‘  *  ,  t  I  Me  (  h  I  ,  ’  I  lOtilou  L  ~  I  IKi'  L  I  ,i  of  the 

i  prii.nv  m  ’  .  I  nor  I. i,i  j  jii  i  tudo  of  the-  tocor.-i ar  y  Mass*’, 

:  in  i  .  .  '  irnqth  c  trio  |  i-  inijry  syr  tom-  '  ,  1  ,  '  the  corf  Me  lent 

:  of  i  -  -  Min  i  on  *  ,  <  .  ■  ' , n ,J  r  y  ■  .  ■ . , . 

:  In  i  i  j  1  ;/■■-'  .v  .1  vo  >  ,  'iho  - 1 .,  i  •  retting*-'  ,  (i  ,  '  Ha-;  I  red 

:  nui.il,  :l  oi  i .  |  1 1. s  -  '  ,  II 

qtv/o  -)j  .  J 

The  I  ni|  •  i  c  t  i  rig  mars  Is  given  an  Initial  position  next  to  the  stop 
oi  poa  I  '  i  l.ir  d  n  vi.tlon  "i  the  Initial  r|.  •  f  1  1.  I  >  n  ei  the  system. 

x  (  0  )  ---t  !  :  t.e  v!  0  )  *1  i  (  d/2  ,  J ) 


Reproduced  Irom 
best  available  copy. 
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* 

* 

* 

* 

* 

* 


pr  1  nt*  , 

•  « 

pr 1 nt* , 

»  • 

print*. 

«  • 

pr 1 nt*  , 

•mass*  * .mass 

print*, 

•  effect  1 vo  d* ’ , d 

de 1 tat«0 

.01 

print*, 

’  e*  '  ,  o 

print*, 

• x 1 0 ) » ' ,x(0) 

print*. 

’daltat*’ .deltat 

pr 1 nt*  , 

•  • 

pr  lot*. 

•  * 

pr 1 nt*  , 

•  * 

The  loop  to  line  900  solves  for  the  motion  of  the  system  through 
n  Impacts. 

do  900  I ■ 1 , n  ,  1 
lnct*=0.01 
da  1  tat  *0 .01 

Here  the  motion  of  the  system  Is  solved  for  after  the  known  time 
of  Impact. 

a( 1-1 >*tdotv( 1-1 )/qone 
b(  1-1  )  «=theta v <  1-1  ) 

thetr.m (  1-1  ) *  sqr t ( a (  1-1  )  *a  (  I - 1 ) +b  <  1-1  )  *b  (  I“1  )) 
aa*(ve1v( 1-1  )  ) 
bb*x  (  1-1  ) 

The  loop  to  line  500  Iterates  to  the  time  of  the  next  Impact  using 
the  requirement  that  the  Impacting  mass  must  remain  between  the  stops. 
An  Impact  Is  considered  to  have  ocurred  whenever  the  Impacting  mass 
comes  within  0.00000 1  feet  of  either  stop. 


do  C00  c=l .500, 1 
time!  1  ) = 1 1  me (  1-1  )+dcltat 

thet  i r I  I  )  *  (  j  (  I  - 1 > "a  I n ( qone*dol ta t ) ♦b { l-l)*cos( qone*del tat ) ) 
x<  I  )  =  I aa  *dc 1  tat )  +bb 
1  1  = theta  s (  1  )  *  1  -  I  d/2 ) 

I  r  ■  1 1 « cJ 

If  (ab'.lxl  I  > -  1 1  >  .UE.  0.000001)  then 
goto  601 

elsclf  (  a  b  s ( x {  I )  —  1 r )  .LE.  0.000001)  then 
aoto  601 

olsc-if  <<x<U-ll>  .LE.  0.0)  then 
lnct-C.5*tnct 
doltat:!doltat-1nct 

goto  560 

elself  llx(l)-lr)  .GE.  0.0)  then 
I  no  t  =  £' .  5  *  I  net 
doltat~doltat-1nct 
goto  500 

else 

doltat=deltat+lnct 
goto  500 
end  I  F 

500  continue 

600  continue 

601  continue 

If  (c  .go.  500)  then 
goto  9'J3 
end  1  f 

If  ( .  L  C  .  2  )  then 
goto  C00 


end  1  f 
contl nue 


The  conditions  Immediately  before  the  time  of  Impact  Just  Iterated 
to  are  solved  for. 

tdots! 1 1  »qone* ! a ! 1-1 )*cos!qone*deltatl-b( 1-1 1  *s  1  n ( qone*del tat  1 ) 
thetav! 1 )*thetas! 1 ) 
ve I s (  1 ) ■ vc  1  v  <  1  -  1  ) 

tdotv (  1 )  =  ( 1 . Jb / (  lnort*nacs*l*l >  1  *  ( (tdots (  1  ) * (  tnert-mass*l*l*ol 
t  +vc  1  s  (  1  )  *mai  w  *  1  *  (  1  .0  +  e  >  >  ) 

velv<  1  1  =  (  1  / <  !nert+mass*l*l > l*(tdots( 1 > 
i  *  I  ncrt*  (  1  .0  +  e  )  + 

i  velst  1  )*ti:idss*l-(  lnert*o/l  )  )  ) 

cont I nue 
goto  994 
cont i nue 

print  ’ < a , 13 1  * , ’ the  time  Iteration  did  not  converge  for  1«',1 
cont 1 nue 

Here  the  motion  of  the  system  Is  given  for  the  first  15  Impacts. 

print  *,' Impact  time  deltat  thetas  x  tdotv  xdot 

:  the  t  un  dthotain’ 

do  9LJ  J-l , 15 . 1 

print  1 ( 17 .  fO . 5,f8.C.f0.4,f9.6,f9.4,f9.4 ,f7.4  ,f9.4>  *  ,  J  , 
i  1 1  iif  -  (  J  <  ,  1 1  J  1  -  t  1  wut  J-l),tliotas(J),x(J>,tdotv(J),velv(J), 

i  th  tii.ilj),tl>otani(J)-thctam<J-l>  . 

con  1  I nuu 
print',  ’  ’ 

print*,  *  ’ 

print*,  '  1 

•  • 

The  errors  in  position  and  velocity  aero ss  the  first  15  Impacts  are 
given  below. 

do  97/T  j»l  ,  15. 1 

err  mom* Inert*! tdotv ( J ) -tdots  <  J) l+mass* 1  * ( vel v( J)-vels(J) ) 
er r  ve 1 re*  <  vo 1 s ( J ) - tdots ( J ) *  1 ) - ( tdotv ( J  >  *  1  - vo 1 v ( J  >  > 
print*.  ’ Imp act*', J. ’ er  r  mom* ’ , err mom, 'errvel*' , err vel 
cont 1 nue 
print*,  ’  * 

print*,  '  * 

print*,  ’  ’ 

print*,  'Impact  thetnm  time  dthetam  dtlme* 

1  th,_  La  *0 . 1 

1 1 1 1110*0 .  0 

de 1 1  j  t =0 . 0 
m*  1 

tt=0.0 
tts-0.0 
mtt-'J.0 
nitts  ;0 . 0 
at*0 . 0 
tat*0.0 
or i  c. r  (01*0.1 
zet;  L “ 0 . 0 
t  z  <  •  v  a  =  0 . 0 
dt(i  )*0.0 
mm*  1 

In  the  loop  to  line  900  the  maximum  positive  amplitude  the  reaches 
during  melt  cycle  Is  obtained  from  the  motion  of  the  system  already 
obtained.  Information  needed  to  do  a  lcjst  squares  to  these 


« 


amplitudes  Is  also  obtained.  The  loop  also  looks  ahead  to  see 
when  the  system  goes  through  a  cycle  without  an  impact  on  the 
assumption  that  thl3  when  the  damper  becomes  Inoperative. 

The  loop  Is  exited  before  the  damper  becomes  Inoperative. 

do  920  JM  ,n-l  .  1 
If  (b(J)  .eq.  0.0)  then 
goto  9U1 
end  I  f 

delt.it=(atanla(J)/b(J)  )  >/qono 
If  (del  tat  .  1  c- .  0.0)  then 
deltats:deltat+(3.14n)<j;:7/qon*} 
end  I  f 

If  ( ( b ( J ) 'cos ( qone*de 1  tat ) )  .)e.  0.0)  then 
goto  979 

elsolf  I  (time!  J  +  12)-ttnic(  J  +  U  >  >  .go.  0.25)  then 
goto  931 
end  I  f 

If  (  ( t  line!  J  )*doltat )  .lo.  time(J+l))  then 
If  I  dil  tat  .go.  0.0)  then 
d  i-.het  i  =  t  hot  unit  J  )  -  1  theta 
t.  t  i  iiio  ( j  )  +do 1  tat 
dL( J )=t-ltlne 

f.  r  i  nt  ’  <  13  .  1‘J  .4  ,f  9  .  4  ,f9.4  ,f  9.4  )  *  ,  J  .thetaml  J  >  ,t, 
i  dtheta , dt ( J ) 

r  es  I  d=*  J 

I  theta9 thetaml  j  ) 

1 t I me»t 
lmpt(m)“t 

r.i  1  nipt  (  mm  )  =*t 

I I  •*  tt  +  t 
t(.3-ttS*t*t 
al»ot ♦thetaml J  ) 

ttat  +  t*thetam(J) 
error  Cm ) 9 thetaml J ) 
n;  =  m+  1 
end  I  f 
end  I  f 
cont I nue 


cont  I  nue 
cont i nue 
m  =  m  - 1 

If  <  ttr.  .eq.  0.0)  then 

goto  9C3 

pis  -If  I m  .eq.  0)  then 

got  o  9 1 3 

els.,  if  I  tts  .eq.  tt*tt)  then 
goto  933 
end  l  f 

A  lci't  squares  approximation  Is  fit  to  the  maximum  amplitudes 
below.  1  ho  rr...  ii.iuiii  d.-v  1 .1.  Ion  from  this  approximation  Is 
obtained,  as  well  as  im.  variance. 

r.«<  I  I  1  .0*m  >  *t.at  )-(  at  *tt )  !/(  I  (  1  .0*m  )  *tts  )  -  I  tt*tt  )  Y 

q  =  t  .  t- I  *lt )  )/(  1  .  0  '  m  ) 

I, ID  .  r  .  0 
v  a i  1-0.0 
do  93.'  J  *  1  ,  m  ,  1 

error ! J  )  *er  r or  I J ) -q - I s" I mp 1 1 j ) ) 

If  labs (error ( J > )  .  ge.  maxorr >  then 
in  l  ■  . ■>  r  ‘irror  I  j  ) 

1 nip  or r  -  j 


112 


ond  1  f 

var I *va p l*orror( J  1  *orror ( j  > 

902  cont  1  nue 

print*,  'up  to  the ' , m , ’ poak  the  least  squares  fit  to  the 
i  peaks  Is  thc-tam*q  +  s"t  with  qo’.q.'and  s»’,s,'wlth  max  error*’, 
t  inaxerr.'at  peak  * ’ , I mpcr r , ’ and  vir i ante* ' , var I 
983  continue 

print*,  ’  ' 

print*,  ’  ’ 

print*.  ’  ’ 

If  I  res  Id  .ge.  n)  thon 
goto  1002 
end  1  f 

If  (resld+25  ,ge.  n)  thon 
goto  1951 
end  I  f 

print*.  '  * 


*  Hero  the  conditions  at  the  Impacts  beginning  where  the  damper  was 

*  assumed  to  bc:o:.w  ineffective  are  given.  This  gives  the  residual 

*  amplitude.  After  that  the  first  10  cycles  of  the  1 nef f oct I vo 1 y  damped 

*  portion  of  the  system  are  given. 

print  *,’ Impact  time  deltat  thetas  x  tdotv  xdot 

s  thrt-.i'i  dthotam* 

do  10 Old  J-r  c-sid. resld  +  25, 1 

print  ’(  I7.f a.0,f8.6,i  1) .  4  ,  f  9 . 6  ,  f  9 . 4  ,  f  9 . 4  ,  f  7 . 4  ,  f  9 . 4  1  ’  ,  J  , 
i  1 1  we ( J ) . 1 1  me ( J ) -t I me( J -  1 ) , thetas  <  J ) , x ( J  > , tdotv ( J  > , vel v ( J ) , 

i  theta  in  (  j  )  ,  thetanil  J  )  -thotuml  J-  1  ) 

1950  continue 

1951  continue 

print*,  '  ' 

print*.  ’Impact  thetam  time  dthotam  dtlme* 

tat*0.0 

lthota-0. 1 

1 1 1 mo*y  .  0 

deltat«0 .0 

ai*  1 

tt-C.0 
tti.*l/.0 
at^0 .0 
mtt=ff .0 
mtt  s  ->0 . 0 
or r o;  <01 *0 . 1 
zet.i  t=0 .  0 
t  zet  a*0 . 0 
dt (0)«0.0 
mm  -  1 

err  <  r ( 0 1 *0 . 1 
do  lilt!'  J»ro5)d,n,l 
If  ( b i J  )  . eg .  0.0)  then 
go  to  19fl 
end  I  f 

deltat“<atan(a<  J)/b<  J) ) ) /gone 

If  id.  I  tat  .lo.  0.0)  thon 
or- It  .  t  -  d«_  1  ta  L+ ( 3 .  141.,527/qone) 
end  1 1 

If  (  (  b  ( .1 )  ’cos  (  qone*dc- 1  tat )  )  .le.  0.0)  thon 
goto  1979 
end  .  f 

If  I  I 1 1  mol J ) +do1 tat )  .lo.  tlmc(J+l))  then 


If  (deltat  .ge.  0.0)  then 

dt  he  t  a  *>  t  ii^tai.il  J  )  -  1  theta 
t  =  1 1  n...-  (  j  l  >  d  .•  1  ta  t 
dt(  J  )-t- It  line 

print  ’<l3.f9.4,f9.4,f9.4,f9.H>’  ,  J,  thetamt  J  ) ,  t, 
t  cl ..  h  c  t  a  ,  j  t  (  J  > 

1 1  h  e  t  a  *•  1 1 1  t  a  in  (  J  ) 

1  t  I  IOC  -  t 

i  nip  t  (  111  )  “t 
ni  I  ript  (  fun  )  ■  t 
tt*tt+t 
tts»tts+t*t 
l  L  ■  1 1  ♦  t  h  e  t  a  m  <  J  > 
tut-tnt*t' uhetam(j) 
error ( m ) *thetam( J ) 
m  -  m  ♦  1 

If  (hi  .  ge  .  16)  then 

goto  1‘JCl 
end  I  f 
Ond  I  f 
end  ( f 

1979  continue 

1 900  cont I nue 

1931  cont I nue 

m*w*  1 

If  <tt«  .eq.  0.0)  then 

goto  1983 

else  If  ( m  . eq .  0)  then 

goto  1903 

else  if  (tts  .eq.  tt*tt)  then 
goto  1983 
endlf 

s=(  (  ( 1 .C*m)*tat)-lat"tt) >/<  < { 1 ,0*m > *tts ) - ( tt*tt >  > 

q»(?t-(s*tt> )/( 1 .0*m) 

maxcrr-W.JU 

var 1=0.0 

do  1333  J* 1 , m , 1 

err or ( J ) “or r or  1 J )-q- <  s*  t mpt( J ) ) 

If  ( abs ( error ( j ) )  .ge.  maxerr  1  then 
max.  r  r  :*er  r  or  (  J  ) 

Imp err « j 
end  I  i 

var  I  *  var  I  *or  r  or  (  J  )  “or  ror  (  J  1 
1  983  cont  I  i.uc 

print*,  ’up  to  the’ ,m, ’peak  the  least  squares  fit  to  the 
s  p.:.il._  I :  tl'iot  ni“c|+o*t  wlui  t)»'  ,i|,  'and  S“’,3,’wlth  mux  error"’, 
i  maxerr  , ’a  L  pc jk -  ’  ,  I inpcr r  ,  ’ and  var  1  ones*  ’  .  var  I 
1933  continue 

999  continue 
1  Oil?.  cont  I  nue 

1000  continue 
1  <30  1  contl  ixie 

end 


Initial  phase  angle*  .000000/T0.Cn+00 

initial  i.. i Ir.ium  deflection*  1 00000001 e+00 

natural  frequency*  .  250000000  --1 02  the  starting  time*  .  0000  100008+00 
moment  of  Inertia  of  the  pr  li.ury  mass*  .  1  09 30000  1  e  +  00 

magnitude  of  the  secondary  mas-s*  .  200..'jLi0009o-02  length  of  the  primary  system 
2.23000002  the  coefficient  of  restitution-*  . 500000O00e+00 

secondary  masses  Initial  velocity* 

. 000000u00e*00  the  gap  setting*  .21100000000+00 
desired  number  of  I nipacts*  50 


mass-  .200000009e-02 
offoctive  d=  . 250000000e+00 
e*  .5000000000+00 
x<0)*  - .9O0000049e-01 
del  tat*  . 999999970e-02 


the  time  Iteration  did  not  converge  for  1*  45 


tmpac  t 

1 1  wo 

de  1  tat 

thetas 

X 

tdotv 

xdot 

thetam 

dthetam 

1 

.06709 

.067007 

.0121 

- .098000 

2.2958 

7 . 0866 

.0926 

-.0074 

2 

. 12207 

.05-130  3 

.0923 

. 330095 

.4363 

-2 . 7636 

.0940 

.0013 

3 

.  10700 

.0665''6 

.0110 

.149614 

-2.2510 

-6 . 2339 

.0907 

-  .0033 

4 

.20345 

.0/5702 

-.0009 

- . 323105 

.2110 

4.0967 

.0093 

-.0014 

5 

.41144 

.  14  790Q 

.0  709 

.2031  1  0 

-1.116  9 

-6.051  1 

.0030 

- .0065 

6 

.50692 

.090100 

-.0031 

- .310375 

-  .  4  4  93 

1 .7270 

.0050 

.0012 

7 

.55043 

.040'.  05 

- .0462 

-.226910 

1 .7234 

4 .9069 

.0025 

.-.0026 

a 

.60214 

. 103713 

.0746 

.291290 

-  .6468 

-4.9173 

.0709 

- .0035 

9 

.70197 

. 119024 

-.0775 

-.297915- 

.  1760 

3.2623 

.0779 

-.0011 

10 

.€0711 

.  005  1  46 

.04  70 

-.020143 

1  .5448 

3.6436 

.0776 

-.0002 

1  1 

. 93917 

.072054 

.049.. 

.235130 

-1  .2663 

-6.2654 

.0708 

•  - .0069 

12 

1  .02/07 

.(10  1503 

-  .  bo  /  6 

-.275461 

- . 7037 

.  9609 

.0751 

.0023 

13 

1  .07296 

.062:.  3S 

-.0440 

-.224004 

1 .3695 

4 . 1 C  10 

.0707 

-.0024 

14 

1 . 1 €756 

. 114606 

.0600 

.264357 

-  .  7969 

-4 .9939 

.0662 

-.004  6 

15 

1 .29300 

. 10551  1 

-.0662 

-.272560 

- . 1705 

2.1134 

.0665 

.0004 

Impact- 

1  errmorn* 

.  ICG? f> 451 5 r* 

-07 

errvol* 

.2394185790-06 

Impact* 

2  orri.iom* 

.  74  6060*, 60o 

-w  J 

e  r  r  v  <_■  1  *  - 

.23H418579e-06 

Impact* 

3  err mow* 

.  7  4  ’.,1 ' 6 0  J  00 ■  ■ 

-III! 

or  r V" 1  *  - 

. 6536743 lCe-06 

Impact* 

4  er riin.im* 

.  0001. 000,. 8.. 

'  00 

OITV.  1“ 

.  47".  137  1  58o-06 

1  nip.,  c  t  * 

5  urrii.uin*  - 

.11  1  760 /i  Jc 

(i  / 

orrv'.-l*  - 

. 23. .4 1 O579e-06 

1  mpu  >.  t  * 

6  crrMOK*  - 

.3  72539 1 

-00 

or  r ve  1  * 

.  2 3 0 4  1  U.#79o-0G 

1  inpac  L  » 

7  er ri.ivM* 

.  ■ .  9 .  ■  i '  3  .  *  2  *,  •  * 

-07 

i  r  r  v  e  1  *  - 

.  6  91.0464  4 lie -0  6 

Impact* 

8  er  rt.i ..  i„=  - 

.  /  4  J  ■)  ’ : 

no 

er  r  Ve 1 =  - 

.  23!: 4  1  05  790-06 

1  mp.'i  cl* 

9  errii'-.m* 

.  jv.'o':90: 

-J0J 

c  r  r  vo  1  * 

.330’ 166790- 0G 

1  mp.i  c  L  * 

10  er  rii  ioiii* 

.  K  l,  17  5  0  l) 

.-07 

er r vo  1  * 

-.11  9  20  92  90.  -06 

Impact* 

11  <-.  rrir.c-iii* 

-  .  2900, '322.: 

,  -07 

err  vo  1  * 

. 2364 1 95790-06 

1  n- pact  * 

12  i  rr  1..0.11 

-.111 756  70  ■ 

-07 

er  r  v«- 1  » 

.  2':i:4  1C579O-06 

Impact* 

13  err.  cm* 

-  .  5  5  i  7  9  3 4 

-  J  Li 

or  r  vi  1  * 

4 1  .:‘..7'h-06 

Imp  i *.  t  - 

14  ..  i  r  i mi* 

-.111  7667'#  i 

-  I*  7 

i :  r r vo  1  * 

-  .  2  1  >’  '•  lio.i /Ou-06 

Impact’ 

15  err  mum* 

.00000'. l  c. 

-•  +  00 

errvol* 

.  4  7  0  3  3  7  1 5  0  e - 0  6 

Impact 

thetan 

1 1  mo 

dthetam 

d  i.  1  mo 

O 

,0'*4r 

.  1296 

- .0060 

.1295 

4 

.1  5  ’  3 

.3063 

-.0647 

.2650 

7 

.  00.:  j 

.6446 

-  .  00  60 

.2592 

10 

.67/0 

.  9039 

-.00.18 

.2594 

1  3 

.0  70,7 

1.1632 

-  .  096  9 

.  2593 

Copy  available  to  Dilt- 
permit  fully  lc(T*fclc  v:p.-du 


IS 

.0644 

1.4230 

-.0063 

.2598 

19 

.0594 

1 .6831 

-  .00 5./ 

.2601 

22 

.0526 

1 .9427 

- .0069 

.2596 

26 

.0410 

2 . 4622 

-.0115 

.5196 

up  to 

the  9 

peak  the  least  squares  fit  to  the  peaks  Is  thetam»q+s*t  with  q* 

.  97664  1 608e-0 1  and  s-  -  .  2296 1 7860e-01  with  max  error-  .  840 1  1  6292e-03  at  peak 
4  and  variance-  . 21 4724355e-05 


Impact 

thetam 

time 

dthetam 

dt  1  mo 

26 

.0410 

2.4622 

- .0590 

2.4622 

28 

.0354 

2.7238 

-.0056 

.2616 

30 

.0298 

2.9863 

-.005  7 

.2625 

32 

.0239 

3.2495 

- .0058 

.2632 

34 

.0180 

3.5136 

- .0059 

.264  1 

36 

.0117 

3.7778 

- .0063 

.2642 

38 

.0056 

4.0299 

- .0061 

.  2521 

39 

.0055 

4.2669 

- . 000 1 

.2370 

40 

.004  4 

4.5122 

-.0011 

.2453 

42 

.0025 

5.2522 

-.0019 

.  7400 

43 

.0024 

5.7657 

-.0001 

.5135 

up  to  the  l 1 

peak  the  least  squares  fit  to  the  peaks  Is  thetam*q+s*t  with  q* 

. 643201 84 3c-01  and  s»  - . 1 23 3497 96e-01  with  max  error*  . 957503542e-02  at  peak 
11  and  variance*  . 3S4233722e-03  . 


program  LAESIM 

THIS  PROGRAM  ATTEMPTS  TO  SIMULATE  THE  CONDITIONS  EXPERIENCED  IN 
THE  LABORATORY  WHILE  FOLLOWING  THE  IMPACT  THROUGH  A  SERIES  OF 
IMPACTS . 

I  n teg er  I , m , c , n , J , p , !mperr,mm,y,u 
parameter  (n=500) 

real  ph I ( 0 :  n  >  ,  thetas  I  1  t  n  )  , thetavl 0 : n  >  , tdotsl  1  :  n  >  ,  tdotv ( Em), 

■ 1 1  me ( 0 : n  )  , ve 1 s ( 1 : n  )  .  ve  1  v  (  it :  n  >  ,  cone  ,  c  two  .a(0:n),b(0:n), 
aa.bb.dtlicta, 1  theta , inct, Itime.t, 

Inert  .mats.  1  ,c,  1  1  ,  dd  ,  gone  ,  q two  .  nit ts  ,  I  mpt  (  1  sn)  , 
ln.pt(  1  :n)  ,  tt,Us,at,tat,error  (iJsn)  ,s,q, 
ma  >:e r  r  ,  v  a  r  I  ,  <J t  (  1  i  r:  )  .  in  I  inp  t  (  1  :  n  >  , 
deltat,  I  1  ,  Ir  ,  x  (  0  :  n  )  ,vel  (  0  :  n  )  ,d,thetam(0tn)  , 
er  r mom  ,  errvel  ,dta,res  Id 
print*,  *  • 

print*,  *  ' 

print*,  *  ' 

The  same  parameters  varied  In  the  laboratory  are  varied  here.  Note 
that  tha  diameter  of  the  Impacting  mats  Is  subracted  from  the  actual 
gap  to  give  an  effective  gap.  Also  note  that  the  same  damping  constant 
Is  not  used  for  each  mass. 

do  1001  p  = 1 , 3 , 1 
If  (  p  .  i  <i .  1  )  then 

mass  =  0 .  P  004  U 1 
ctwo-U  .  0/1.000  I  7 
d  la-7.  (1611/12.0 
el se if  (p  .  eq  .  2)  then 
mass  --0 . 001  4  9 
ctwo=0 . 011002  35 
419*0.6875/12.0 

else 

masr*0. 00503 
c  two=0 . 00002  67 
dla-I .031/12.0 
end  I  f 

d*  3 .0*  < 1 .0/12.0) 
d“d-d  I  a 
e*0. 5 

Previous  values  of  the  system  position’s  and  velocities  are  set  to  zero* 
do  1C1  J* 1 , n , 1 
thetas! J  >=0.0 
thotavl j )=0.0 
tdotsl J ) =0.0 
tdot v ( J ) =0 . 0 
time! J ) =0.0 
vel s<  J  >*0 .0 
ve  1  v  (  J  )  =0 .0 
x<J)=0.0 
thotaml 3 ) *0 .0 
cont I  nun 
print*,  ’  ’ 
print*,  '  ’ 

print*.  ’  ’ 

The  remainder  of  the  system's  parameters  are  defined. 

data  ph 1 ( 0  )  .  thetav (0 ) .t lme<0) , qone , Inert . 1 . 

>  vc  1  v  (  0 ),  in  ,  tdotv  ( 0  >,  cone  ,  1  1  / 


I 


0.0,1 .0.0,0.02,7 .0/ 

pr Int* , ’ In 1 1 1 a  1  phase  angle"',  ph 1 10) , 'Initial 
maximum  def 1 oct 1 on= ' ,  thetav(O), 

'the  starting  1 1  me* ' , 1 1  mo < 0 ) , ' moment  of  Inertia  of  the 
primary  mas s * ' , I ner t , ’ inagn I tude  of  the  secondary  mass"', 
mass, 'length  of  the  primary  system" ', 1 ,’ the  coefficient 
of  restitution*’ ,e, 'secondary  masses 

Initial  veloc 1 ty" ' , wel (U ) , ' the  gap  sett  I ng« ’ ,d , ’ des I  red 
number  of  Impacts"’ ,n 
qono*25 . 9 

qtwo  =  sqr  t (<32.l74/ll)-( c two/ ( 2 . 0*ma s s  ) >*( ctwo/ (2.0*mass))> 
print*,  ’  ’ 

print*.  '  ’ 

print*,  '  ' 

The  Initial  position  of  the  Impacting  mass  Is  placed  against  the  stop 
opposite  to  the  Initial  deflection. 

x  (  0  )  “theta  v(0)*l*(d/2.0> 
print*,  ' ma s  s ■ ’ , ma s  s 
print*,  'effective  d"',d 
del  tat  =£T.ill 
print*,  ’ x ( 0 ) ■ ’ , x ( 0) 
print*,  'deltat* ' .doltat 
print*,  ’©■ * ,o 
print*,  'gone*', gone 
print*,  '  qtwo"  ’  ,  qtwo 

The  loop  to  line  9£(0  solves  for  the  series  of  Impacts. 

do  900  1  *  1 , n , 1 
I nc t-0 .01 
dcltat=0.01 

The  motion  Is  solved  for  after  a  known  time  of  Impact. 

a ( 1-1 ) ■ ( tdotv ( 1-1 )/qone+thetav< 1-1 )*cone/<2.0*1 nor t*qone ) ) 
b<  1-1  IHhc-tavl  1-1  ) 

thettral  i-1  )=sqrt(a(  1-1  )*a(  1-1  )*b(  1-1  >*b(  1-1  )  ) 
aa *  I  v<_  1  v  l  1-1  ) /qtwo+ ( ctvo'xl  1-1 ) / ( 2 .  0*  lira  s  s  *q  two  > ) ) 
bb  =  x (  1  -  1  > 

The  loop  to  line  GOO  Iterates  to  the  noxt  time  of  Impact.  This 
Iteration  uses  the  requ  1  r cuc-nt  that  the  Impacting  mass  remain 
between  th-  stops.  An  Impact  Is  d'.  f  1  nod  as  ocurrlng  whenever 
tlio  Imparting  mass  Is  within  0.U0OOOI  foot  of  either  stop. 

do  600  c"  1  , 500, l 
t  Im*  (  I  l^tl  ll»t_  (  I  -  1  )  +dc 1  tat 

t  he  t  ,s  (  I  )  =  (._>:  p  (  <  -(  cone  )/  <  2 . 0*  1  ner  t )  *do  1  tat  )  )  )* 

:  <-(1-l)*sln(qC'rio*deltat)+b(l-l>*i:o5<qone*deltat)) 

x  (  I  >  =  (cxp(  (  -  (  ctwo  )/(2.u*mjaJ,))*doltJt)l* 
s  ( a a  *  s 1 n ( q two* do ltat)+bb*cos( qtwo* do  1  tat ) 1 

11»thetas(  I  )*l-(d/2) 

I  r  »  I  1  *  d 

If  ( abs I x I  11-11)  .LG.  0.000001)  then 
goto  601 

elsolf  ( abs ( x (  1  )  -  1  r  )  .LE.  0.000001)  then 
goto  60 l 

el  so  if  ((x(l)-ll)  .LE.  0.0)  then 
<  n c  0 . 5*  1  net 
do  I  t  :i  I  •  de  1  tat-  1  net 
goto  600 


elself  ((x(l)-tr)  .G£.  0.0)  then 
1 nct»0 . 5* I  net 
de)tat»deltat-lnct 
goto  500 
else 

de)tat*deltat+1nct 
goto  500 
end  1  f 
cont  t  nue 
cont I nue 
cont 1 nue 

If  (c  .gft.  500)  then 
goto  993 
end  I  f 

If  (c  .IE.  2)  then 
goto  800 


2)  then 


end  I  f 
cont I nuo 

The  conditions  Immediately  before  the  time  of  Impact  Just  Iterated 
to  are  solved. 

tdots ( I )“-(cone/(2.0*inort) >  * (exp ( (-cone / ( 2 . 0*1 ner t >  *de 1  tat ) ) ) * 
i  (  r.  (  1  -  1  )  *s  I  n(  gone*  do  1  tat )  +  b  <  1-1  )"cos(qono*deltat))+ 

i  (  e:-  p  (  (  -cone/  >2.0*  lnort)*deltat  )  )  )  "qone* 

i  (a(1-l)*cos(qone*deltat)-b(1-l),'s1n(qone*deltat)) 

i 

thetav<  I  )*thetas( 1 ) 

vel s ( l >■- ( ctwo/ »  2 . 0*mass ) )*(exp( ( -ctwo/ (2.0*mass)*deltat) )  J  * 

«  (aa'j  ln(  qtv/o*de  1  tat )  +  bb  *cos  (  qtwo*de  1  tat )  )  + 

i  ( exp ( ( -c  two/ <2.0*mas3)*deltat) >) *qtwo* 

i  ( a a  *cos ( q two 'do ltat)-bb*s I n ( qtwo*de 1  tat > ) 

tdotv( 1 )•< 1 .  0/  <  1  nor  t  +  ma  s  s  *  1  *  1 > )*( ( tdots ( I >*( 1nert-mass*l*l*e) 

:  +vel  (  l  )  ‘mass*  1  *  U  .0+e  )  )  ) 

ve ?  v< I )  *  < 1 / <  Inert  +mass*  J  *  I  )  }  *< tdots (  I ) 
i  *  I nert* ( 1 . 0  +  e  >  + 

:  ve 1 s ( 1 ) * ( mass" 1  - < I nert*e/ 1 ) ) ) 

continue 
goto  994 
cont I nue 

print  ‘ ( a , 1 3 ) ’ , ’ the  time  Iteration  did  not  converge  for  1«’,1 

continue 

cont I nue 

Here  the  system’s  condition  at  the  first  15  Impacts  Is  given. 


print  *,’ Impact  time  doltat  thetas  x  tdotv  xdot 

:  Diet  am  <Jt  h..- L.am  ’ 

do  ‘jtl  J=  1  .  1  v  ,  1 

print  ’(  '7,ff!.5,f8.r.  ,fC.4,f9.6,f9.4,f9.4,f7.4.f9.4),,J, 

:  t  I  .,!••(  J  )  .  L  1i,i.  .1  )  -  1 1i>;..-(  j  -  1  )  .  thetas  (  J  )  ,  x  (  j  )  ,  tdotv  (  J  )  ,  vo  1  v  (  J  )  , 

!  t h t in  ■  J  >  ,  1 1 : i- 1  a iii (  J  )  tl)Ct.mi(  J  -  1  ) 

cont  I nuo 
cont  f nue 
print*.  '  ’ 

print*.  ’  ’ 

print*,  '  ' 

The  errors  In  the  momentum  and  velocity  across  the  first  1  Impact: 
is  given.  Ylioy  ..Mould  to  si.iill  enough  to  be  a5:.umod  ne-jllblo. 

do  97U  J  =  l,lr..t 

err mom = I  ner  t  *  ( tdot  v  l  J  )  -  tdots  (  J  )  )  +m.a-.;s  *  1  M  ve  1  v  (  J  )  -  ve  1  s  (  J  )  ) 


tdotv 


I 


119 


9 


errvel =  e*< vclsl  j >-tdots<  J >*1 ) - ( tdot v ( J ) *  1  - ve 1 v ( J  >  ) 
print*,  ’ Imp act*’, J,’ err mom* ’ , er rmom , ’ err ve 1 “  ’ , errvel 
cont Inuo 
print*.  ’  ’ 

print*.  *  * 

print*.  *  * 

print*.  'Impact  thetam  tlmo  dthetnm  dtlme' 

I  theta -0 . 1 

It i mo  =  0 . 0 
de 1 t at *0 . 0 
m»  1 

tt*0.0 
tts=0.0 
at*0 .0 
tat*  0 . 0 
mtt*0 .0 
mtt s  *0 . 0 
c-r  r or  (0)  *0. 1 
dt ( 0 ) *0 .0 
mm*  1 

The  Tol lowing  loop  takes  tho  series  of  positions  and  solves  for 
the  pc?!'  positive  amplitudes.  It  also  looks  ahead  to  see  when  the 
system  rioo;  through  a  complete  cycle  without  an  Impact  to  on  the  that 
this  is  when  the  damper  becomes  I  nc-ffeev  I  ve .  Information  needed 
to  perform  a  least  sviu.iros  fit  to  these  peaks  Is  also  obtained. 

do  930  J*2 , n-2 , 1 
If  (  a  (  ,i  )  .  eg  .  0.0)  then 
goto  900 
end  I  f 

plil  (  .1  ) *  a  ta  n I b (  J  )  /  a  !  J)  ) 

deltat*(atan(2.0*1 nor t*qone/ cone )-phl(J))/qone 

If  (doltat  .le.  0.0)  then 
d«,ltJt*doltat  +  (3.141'i927/qone> 
end  I  f 

If  ( ( b ( J ) *cos ( qono*de 1  tat ) )  ,1e.  0.0)  then 

goto  973 

elaolf  ( ( 1 1  me ( j  +  2 ) - 1 I  me ( J  +  1 ) )  .ge.  0.25)  then 
•  lot o  981 
end  I  f 

If  ( <ttmo( J)+dcltnt)  .le.  tlme(j+i>)  then 
If  (de-ltnt  .mo.  0.0)  then 

bth  .  L  . . .  J  )  -  ltheta 

t  *  t  i  r,  •.  (  .1  )  nJ'.  1  tat 
•J  t  (  J  )  *  t  -  1 1.  1 1,  e 

print  ’<  i  3  .  i  5 . 4  .  f  9  . -I .  f  9 . 4  ,  f  9 . 4  )  ’  .  j  .  thctaml  J  )  ,  t , 
t  d  theta . dt ( J  ) 

r  r.  i  d  *  )  *  2 
lthoti*thotam(J) 

1  t  I  lilv  t 

I  i.ip  L  ( in  )  *  t 

in  imp  t  <  m.ii )  *  t 

tl»tl*t 

i  iti*t*t 

•  it-  l  HI:  t  1.71  <  J  ) 

t  .  I  i  it  1*1  h.-l'iini  j  ) 

<  r  Oi  (  I.. !  *  t  lie  t  am<  J  ) 
ii.  -  in  +  1 
•  end  I f 
en  .If 
cein  1 1  nue 


900 

981 


cont 1 nue 
cont  t  nue 


m*m-l 

if  <tts  .eq.  0.0)  then 
goto  'j? j 

els o If  (m  .eq.  0)  then 
goto  903 

elseif  (tts  .oq.  tt*tt)  then 
goto  903 
end  1  f 

*  The  following  equations  solve  for  the  least  squares  fit  to 

*  maximum  amplitude  of  the  cycles  just  obtained.  The  maximum 

*  departure  from  this  least  squares  approximation  and  the  variance 

*  Is  also  obtained. 

s»( ( ( 1 .0*m>*tat)-{at*tt) )/( < ( 1 ,0*m)*tts >-(tt*tt> ) 
q=(  ;.t-  C  s*tt )  )/<  1  . 0  *  m  ! 
ma  x  - r  r  ~U . 0 
var  1  ~SS  .  £) 

do  S 82  J*1 ,m, 1 

error<j>»error(J)-q-{s*1mpt(J)) 

If  ( abs^ error ( J ) )  .ge.  maxerr)  then 
ma xor r “error ( J  ) 

Hiporr**  J 
end  1  f 

var l=var 1+error( J I'error ( J) 

982  cont 1 nue  ’ 

print*,  'up  to  the m peak  the  least  squares  fit  to  the 
»  peaks  Is  thetam=q+s*t  with  q=’,q,’and  s=',s,’w1th  max  error*' , 
i  maxerr,  'at  peak  * ’ , I mper r , *  and  var lance* ', var 1 
print*,  ’  ' 

print*,  '  ’ 


pr  1  nt*  ,  ’  ’ 

2905  continue 

903  continue 


If  (rcold  .go.  n)  then 
goto  10U2 

e/  if 

It  i  o".  ld  +  25  .ge.  ri )  then 

goL.  1951 
end  1  f 

print*,  ’  ’ 

* 

*  Her?  the  rondltions  at.  tho  Impacts  boqlnnlnq  where  the  damper  was 

*  3’.  d  to  hr... M..:  1  no  I  f  •  ••:  f  I  ve  are  given.  Inis  lives  the  ro.lduil 

ar.plitu  After  th  ,t  the  first  15  cycler,  of  the  I  nof  1  net  I  ve  1  y  d-mped 

*  portion  of  tho  system  ;.re  given. 

* 

print  *, ’Impact  time  deltat  thct3S  x  tdotv  xdot 

:  tlioi.ni.:  *J  1 1 .  .  t  ti  :a  * 

do  1950  J »r es 1 d . r es 1 d+ 25 , 1 

print  ’  (  17  ,f8.5,f8.6,f8.4,f9.6,f9.*,f‘  4  ,  f 7 . 4  ,  f 9 . 4  )  ’ , J , 
i  t f me ( J > , t ) me ( J ) -t Ime 1 J- 1 ) , thetas  1 J > ,  1  > , tdotv ( J ) , vel v l J ) , 

i  thetam(J>,thetam<J)-thetam(J-l ) 

1950  continue 

1951  continue 

pr 1 nt*  ,  *  ’ 

print*.  ’Impact  thetam  time  dthetam  dtfme* 


■  t « 


tat -0.0 

1 thota'0 . 1 
1 1 1  ni.i  =  0 .  0 
del tat=0.0 
m«  1 

tt=0.0 
tt  s  =  0  .  0 
at=0.0 
mtt“0 . 0 
mtt s  *0 . 0 
er  r or  <  0 ) =0 . 1 
z  etat “0 . 0 
tzeta*0 .0 
dt<0>-0.0 
mm«  1 

er  ror  <  0 ) “0 . 1 
do  1980  J»resld,n,l 
If  (a(J)  .eq.  0.0)  then 
goto  1980 
end  1  f 

phl(J)=atan(b(J)/a(J)) 

d<iltat=(atan(2.0*tn<-.rt*qone/cone)-pht(J))/qone 

If  (dcltat  .le.  0.0)  then 
do  1 1  a t*de  1 1  a  i.*  (  3 . 1  4  1  092 7/  qone  ) 
end  1  f 

1 1  ( ( b <  J ) *  cos ( qone*dc 1  tat ) )  .lo.  0.0)  then 

goto  1979 
end  I  f 

If  ( ( t I me( J > +oel tat )  .le.  tlmelj+l))  then 
(f  (dcltnt  .go.  0.0)  then 
dthct  a  =  lliOt  i in (  J  )  -  )  theta 
t  ■  t  i  i...j  (  j  l  *  de  1  tat 
o  w  l  j  )  - 1  -  1 1 1  ni.  - 

print  ’(  I3.f9.4,f9.4,f9.4,f9.4)'  , J.thetaml J )  ,t . 
t  dthcta.dt(J) 

1  thet  a  *  thetani  (  j  ) 

)  1 1  c.ie»t 

I  i.ip  l  ( in  )  =  t 

(,i  1 1 1;  ■  L  ( i.iiii ) 1 1 

I I  -  1 1  ♦  t 

It  ;  =  tts.  +  t*t 
a  t 1  ♦  1 1 1  <  t  a  n  (  J  ) 
tat  t.it  H  Hlmlaml  J  ) 
t.  r  r  n  r  (  in  )  ••  1 1 1  e  tai.il  J  ) 

in  "  hi  ’  1 

II  (in  .  n  e  .  18)  then 

goto  1981 
end  I  f 
end  I  f 
end  I  f 

1979  continue 

1  ‘..80  coil'  I  nue 

1 9 u 1  con  t I nue 

ip'r.i-  1 

If  (tt:;  .eq.  0.0)  then 

not..  1983 

cl.'.iV  (m  .cq.  0)  then 

g>»i.o  1983 

el  •  If  (Li  .  .eq.  tt*tt)  then 
goto  1983 
end  I  f 

oM  i  (  1  .C*ni>*tat  )-(nt*tt)  )/(  (  (  1 . 0*m )  *  tt  s  )  -  ( tt*  tt )  ) 
q*(Ll-(5*ttl )/( 1 .0*m) 


I 

I 


s 


maxcrr  *0 . 0 
var 1-0.0 
do  l‘Jd2  J* 1 , m . 1 

err  or  <  J ) 'error ( J ) -q- ( s* I mpt ( J )  ) 

If  (abs,  (  error  (  j  )  >  .go.  maxerr  )  then 
inaxerr “error  (  j  > 

I  taper  r  =  J 
end  1  f 

var  l=varl+orror(J)*error(J) 

1982  continue 

print*,  ’up  to  the' am, ’peak  tho  least  squares  fit  to  the 
i  peaks  is  th;tai..- g  +  s*  t  with  q=’,q,'and  3=* '  ,  s  .  ’  w  I  tli  max  error*' 
t  maxerr, 'at  peak  =  *  ,  I  mper  r  ,  ’  and  var  1  anr.e=  '  ,  var  I 

1983  continue 
999  cont I nue 
1002  continue 
1 IJ00  cont  I  nuo 
1001  continue 

end 


Initial  phase  angle®  .13080000000*0111 

initial  maximum  deflection®  - . \8888880l e*80 

the  starting  time®  . 000000000e+00 


monc'it  of  Inertia  of  the 

magr.  I  tude  of  the  secondary  mass® 

2. 2 1000004  the  coefficient 
so' ondary  masses 
. 01007C000C < fie  the  gap  setting® 
desired  number  of  Impacts 


primary  mass®  .  1  G7999994e  +  00 
.  4009999  'JCe-i  .i  length  of  the  primary  system* 
of  restitution-  .  5h'0o08000e+00 
Initial  velocity® 

.  21091G668ef)JU 


mass®  . 4O0939995e-03 
effective  d®  .  2 109 1 G668e+00 
x(0>®  -.11094101.70  +  03 
del  tat®  .9939999700-02 
e®  .  00.900'U +  00 
qonc*  .  200 9 9 9 s9Ge+02 
q  two  ®  2.  14377C42 

the  time  iteration  did  not  converge  for  I1 


ct 

t  1  no 

d  c  1 1  a  t 

thetas 

X 

tdotv 

xdot 

thetam 

dthetam 

1 

.05912 

.059122 

-  .00  4  1 

- . 1 1 4G1  G 

2.5324 

0 . 431G 

.0979 

-.0021 

o 

i. 

. 10972 

.050594 

.05  32 

.31 1323 

.  3051 

-1  .5881 

.093  2. 

.0004 

3 

.  158.9  1 

.  04  In  3  5 

.  0  5  7  2 

. 23 1 0G5 

-2.0393 

-5.9662 

.0972 

-.0010 

4 

.  2  50  5  9 

.091042 

-  .0952 

-.315095 

.  3935 

4 .3304 

.8904 

- .0008 

6 

. 308  1  7 

.13  /  5  1  3 

.0002 

.282787 

-1  .  2935 

-6.51 24 

.0944 

•- .0020 

6 

.  4  7.. -.8 

.091:>13 

-.0912 

- . 30708 1 

-.6221 

1  .  22G  1 

.  0944 

.0000 

7 

.51/42 

.113  8  7  4  1 

-.0093 

-.250617 

1 .6380 

4 . 8002 

.  0937 

- .0007 

8 

.51/3  12 

•035703 

.0080 

- .006497 

2 . 4007 

5.5004 

.0933 

- .0003 

9 

.  G  2  8  7  7 

.009041 

.0081 

.300210 

-.7131 

-5. 1092 

.0923 

-.0011 

10 

.73913 

. 1 10302 

- .0503 

- .305007 

.  3b3  1 

3 . 8401 

.0914 

- .0009 

11 

. 80337 

•0G423C 

.0221 

-  .050052 

2.2773 

5.6191 

.0907 

-.0007 

12 

. 8GG31 

.002947 

.0802 

.295993 

-.6490 

-5.0017 

.0097 

- .0010 

13 

.98481 

. 1 18496 

- .0072 

-.290257 

.4319 

3.9743 

.0008 

- .0009 

14 

1.12431 

. 139504 

.0097 

. 259508 

-1.3413 

-G.4926 

.0068 

-.0020 

15 

1.20832 

.084007 

-.0018 

- . 286139 

-.7575 

.  7055 

.0009 

.0001 

Impact® 

impact* 

I mpa  c  t  * 

Impact® 

I mpa  c t  = 

Impact® 

Impact® 

Impact* 

I mpa  c  t * 

1 mpac t ® 

1  mpa  c  t  = 

I  mpa  c  t 

Impact® 

Impact* 

Impact® 


1  errmom*  - 

2  or t mom*  - 

3  errmom® 

4  errmom®  - 

5  cr  i  i  .-'in* 

G  errmom*  - 

7  errmom*  ■ 

8  errmom*  - 

9  errmom* 

10  or  r  I'/.oiii* 

1 1  e r r mom* 

12  err  ......I * 

13  cithioiii* 

14  err 

10  err  ii.oin* 


4004 6 8 70 7 e -07 
2  32C  201.4  <1  e-07 
2923uG01 1 e-07 
3725290300-03 
14  901  It.  i  ..‘e-07 
331  32257 5o-09 

20*34203::.!.;  -07 
55  1.3010000-0  7 
22301 74 1 8e-07 
.C:.  C  19;J2  l  7. .--00 
.4/1  4 o:.*  32... -07 
.9  1 322!  7ho-09 
. 1 ..02045 1  he -90 
.  9  2  1  32  2  0 7 Oe-0<J 
. 1 30921 128e-07 


errvel-  . 1 1 9209290e-05 
errvel®  . 7 1 5255737e-06 
errvel®  - . 71 5255737e-0G 
errvel*  . 4 7 0 0 3  /  1 5Ge-0G 
errve-1®  .  0000 00000 e  + 00 
errvel*  - . 2384 1 OO79e-06 
errvel®  . 2301 18G79o-0G 
errvel®  . 9 03b 7 4 3 1 6c-0G 
errvel®  -.47bU371O3e-06 
errvel®  .  00P0ti0000e  +  00 
errvel®  .  9  0  2.  o  7  I  'a  1  G  e  -  0  6 
err vel »  - .238 4 IU5/9C-0G 
errvel*  . 2304 1 O579o-06 
errvel®  - . 4 7683 7 1 5Ue-0b 
c.rrvel®  - . 7 1 5255/ 37e-0G 


Impact  thetam  time  dthetam  dtlme 
2  .0983  .1221  -.0017  .1221 

4  .0101  . 3  G 1. 8  -.0019  .2437 

8  .0933  .0101  -.0031  .2443 


penult  k»  ,Gl  =  P 


124 


1 1 

.0907 

.0544 

-  .002G 

.2443 

13 

.0883 

1  .  0  9  „  8 

-.0019 

.  2444 

1  7 

.0868 

1  .  3431 

-  .0029 

.  24  4  3 

20 

.0834 

1  .  56  7  6 

-  .  Oj-25 

.2443 

23 

.0811 

1.8319 

- . 0023 

.  2444 

26 

.0703 

2.0763 

-  .0023 

.2445 

23 

.0769 

2 . 3206 

-  .0019 

.2444 

31 

.0744 

2.5651 

-  .0025 

.2444 

34 

.0720 

2  .  U0  9  6 

-  .6024 

.2445 

36 

.0702 

3.0541 

-.001*3 

.2445 

39 

.0677 

3 .2986 

- . 0025 

.  2444 

41 

.0660 

3.5431 

-  .  OX*  1  7 

.2445 

44 

.0635 

3.7875 

- .0025 

.2444 

46 

.  OG  1  9 

4.0321 

- .00 1 6 

.2446 

49 

.0593 

4.2763 

- . 0O26 

.2443 

51 

.0570 

4 . 5210 

-  .01*16 

.2446 

54 

.0553 

4.7653 

-  .0025 

.2443 

56 

.0539 

5.0100 

-.0U15 

.2447 

59 

.0515 

5.2542 

- .0024 

.2442 

62 

.0495 

5.4907 

-.0020 

.2446 

65 

.0477 

5.7434 

-.0010 

.  2446 

6Q 

.0450 

5 . 9070 

- .001 0 

.2445 

72 

.0422 

6.4766 

-  . J'l*30 

.4008 

74 

.0  406 

6.7213 

-.1*018 

.2446 

76 

.0386 

6.9659 

-  .0017 

.2447 

78 

.03  7  1 

7 .21x17 

-  .00 1 7 

.  24  48 

lib 

.0354 

7 . 4555 

-.0017 

.2448 

02 

.0337 

7 . 7 00  4 

-  .00  1 7 

.2449 

84 

.0321 

7 . 94  5  4 

-.0016 

.  2450 

UG 

.08  0  a 

C.  190  a 

-  .  00  1  6 

.2451 

88 

.0280 

0.4366 

-.001  6 

.2461 

90 

.0273 

8 . 60X83 

-  .10  16 

.2452 

ti  *> 

.  0  3  5  6 

0.9261 

-  .Ox*  16 

.2462 

94 

.  H  2  4  2 

9.1714 

- . ro  1  3 

.  24  63 

96 

.0227 

*• .  4  1  G  7 

-  .0015 

.24  64 

98 

.8212 

9 . 6621 

-  .00  1  6 

.2454 

l£C 

.0196 

9 . 9076 

-  .0016 

.  2-155 

10  2 

.  v  1  6  1 

lx.  1  3  : 1 

-.0015 

.24  05 

104 

.0165 

1  0 . 3  9  0  0 

-  .00  l  6 

.2455 

IOC 

.615  9 

1  0 . 644  1 

-  .0016 

.2455 

108 

.0184 

1-'  .  8096 

-  .  0 .  •  l  6 

.2454 

1  10 

.<•.118 

1  I  .  1 346 

-.0xit6 

.  245  1 

l  12 

.  0  1 0  1 

11.3791 

-.001 6 

.24  45 

1  1  4 

.  C  >.t«36 

1  1  .6223 

-  .Ox*  16 

.  2432 

up  to  1.1. u  4  7 

pill  tin  pqu.ir"!  fit  to  the  peak?  Is  thetam=q  +  s*t  with  q  = 

.  J. 5o ; .  325<_  -o  1  u  1 1  >J  v.*-  -  . /O',  x)  1  668  le-02  with  max  error*  , 4/O635900C-02  at  peak* 
2  jrid  v.irliincu-  .  2  x  6  6  G  4  6  7  /  u-  -  0  3 


Impact  t  1  r  -  d  <?  1 1.  a  t 

thot  a  J 

y. 

tdot  v 

xdot 

thetam 

dthotam 

1  l  6  1  1  .  0  3  6  :  3  .  1  3  4  1  1. 1 

.  0)  4  1 

-  .096373 

.  1562 

1  .  7656 

.0073 

-  .01*0  6 

11711.971/4  .14(09 

-.21  X .  9 

.090233 

- .0407 

- . 797  1 

.0071 

- .0002 

1  IC12.2J  .134  .26:069 

-  .  ,0o OU 

-  .  1  7.0 500 

.03/6 

.5171 

.  00  7  0 

- .0002 

119  1.:.  i  :  .411101 

.0X6J9 

.  1.97523 

-  .  16/2 

-.8221 

.  00. .6 

-  .  0004 

1/0  1/.  8  1.0  ...68ir,7 

-  .  *  ‘  •  1 1  7 

- . 100283 

-.1611 

-  .  24  4 

.  /m. 6 

.0001 

1211::.  9/1  ^  .  x*  725  69 

-  ,  .  1 1  0  5 

-.117  6 4 1 

.0898 

.  14  6  9 

.  00  6  5 

-  .  Ulllj  1 

) '21/. 5 V 1 37  .597509 

.  1.1  (.0 

. 1 1  f 1 0  37 

-  .048 1. 

.351/ 

.  00  6  2 

- .0102 

1  »  1  4  . IO/14  .6145u9 

-  . 1".  . 3 

-  .  I  1  /  14X1 

.072  2. 

.•U  4.1 

.60  00 

- .0003 

l.'.i  1  4  .  6  6208  .464  1  4  1 

-  .  t'X'66 

.092606 

- . 0052 

-.2529 

.  00  6  3 

-  .XTO02 

126 16. 41'.  .707514 

-  .15'  J  1 

-.11  2/6  1 

.  1  1  72 

.  4  5P  2 

.0056 

- .0004 

'2615.82'.  <0  . 40JC/3 

-  ,  .0  0  JxJ 

.096901 

-  .  1090 

-  .  Co  a  3 

.0062 

- .0003 

•«> 

• 

* 

• 

12716. 159 

37  , 

,  33,9205 

- . 00 1 3 

- . 108324 

.1195 

.  7005 

.004  8 

-  .0004 

12016.470 

08  . 

.31  11.09 

.  00  4  7 

.  1  1  0304 

.012  4 

-.3870 

,604/ 

-.0001 

12917.106 

.38  , 

,630002 

-.0i,8C 

-.11301! 

.  0  6  7  9 

.3783 

.  004  0 

- .0002 

130 17.61  CIO  . 

.  5  1 1 U  9  0 

-  .  Dif  i  3 

.  102070 

.  lir  Jl 

.lu'.'  1 

.0044 

-  .  0110  1 

13117.69031  , 

,002335 

.in}  4  2 

.  1  1  4  Li  2  U 

-  .  U2U0 

-  .  1  5U3 

.  00  4  4 

.0800 

13217.74047 

.04  1660 

.001  1 

.107/91 

-.  1092 

-.2729 

.  8043 

.01100 

13318.37001 

.638134 

.  0  u  2  3 

-.  100  4 1 0 

.0872 

.  4307 

.  004  1 

-  .0003 

13418.01302 

.430210 

-  .06 ' 2  6 

.1199300 

.0810 

.0050 

.00  40 

.0000 

13010.973 

99 

.  1601  CC 

-.0.0  i  3 

.  102670 

-  .0967 

-.3126 

.  8039 

-.0001 

13619.560 

27 

.  59 1 286 

-  .0./..I3 

- . 106U37 

.  0  9  4  0 

.  4  b  b  1 

.0037 

- .0803 

13719.970 

09 

.  4  106  20 

-  .0632 

.09386  1 

- .0340 

-.3499 

.0035 

-.0001 

13020. 5 1903 

.543130 

-  .  ij  u  i  fi 

- . 100/74 

.0701 

.4108 

.0063 

- .0002 

13920.97719 

.  400 164 

-  .0629 

.098359 

.034  f 

-.0960 

.0032 

.0000 

14022.230 

30  1 

.250192 

.0001 

- . 100289 

.0767 

.  2575 

.00  30 

- .0003 

14122.893 

06 

.  650272 

-  .0028 

.099267 

-.0124 

-  .  1724 

.0028 

-.0001 

Impact  thetam  time  dthetam 

dt  1  me 

116 

.0073 

1  1 

.  8626 

-  .0927 

1  1  .  8626 

118 

.0070 

12 

.  34  70 

-  .0003 

.4044 

121 

.0060 

13 

.0708 

- .0005 

.7287 

123 

.00  60 

14 

.  29. j0 

-.1-1*1  0  sJ 

1.2147 

125 

.0095 

1  0 

.  GjU 30 

-  .  00/*  0 

1.2126 

127 

.0048 

16 

.2304 

-  .0007 

.  7274 

128 

.004  7 

16 

.  4746 

-  .0001 

.  2444 

129 

.0045 

17 

.2036 

-  .000  2 

.  7288 

130 

.6044 

17 

.  fcj  Cl  I  j  G 

- .0001 

.4851 

133 

.0041 

10 

.4162 

-  .0003 

.  7276 

134 

.  0040 

in 

.  90ii9 

.0(9.9 

.  4C47 

136 

.0037 

1  9 

.  6266 

-  .0004 

.7277 

138 

.0933 

20 

.5981 

- .0804 

.9694 

* 

139 

.0032 

e  » 
4-  4 

.0080 

.006  0 

.4839 

140 

.0030 

22 

.2900 

- . 00t  3 

1.2129 

up 

to  the  15 

peak  the  least 

squares 

fit  to 

the 

peaks 

Is  thetam 

»q  +  s*t 

with  q* 

. 

1/0033 150c 

1-01 

and  s“ 

- . 4 2  3  3939  9OC-03 

with  max 

error* 

398508273e-03  at 

10  ani  vartuncc'  .6l7C17292e-06 


Initial  phase  angle*  . C00000000O+00 

initial  Imuni  deflection*  -.1.000000010+00 

th j  starting  time*  . £000000 6 8e *00 

moment  of  Inertia  of  the-  primary  mass*  .  1  O79999‘J4e+00 

magnitude  of  the*  secondary  mass*  .  1 4 -Ji  0000OO-02  length  of  tlio  primary  system* 
2 . 2 10O0D04  the  coefficient  of  rest  I  tut  Ion*  .  O’j000'0000e+00 

sc.'  ■  .n.j  .  r  y  ma ..  r,<.  c  Irilflal  vploclty* 

.01. JWu0000e+tD  the  gap  setting*  .1927003260+00 
deslreu  number  of  impacts*  000 


mass*  .  1  4  Of;0f.:005  e -02 
effective  'l~  .  1  9 2 70 l]328e+00 

x ( 0 ) ■  - . 1 2  I « •  •!  0  8  3 / e  +00 
del  tat*  .  9'  J'."J  9  9  7  6e -02 
c:«  .  0000001  '0  00 

gone*  .  25U9'3'3 aOOe+02 
q  two*  2.14"  >ugi  L  1 

the  time  Iteration  did  not  converge  for  I*  71 
Impact  time  do  I  tat  thetas  x  tdotv  xdot 

1  .001.2!  .  0!ir. n  1 0  -.012  1  -.12270'’  2.4205  C.1C53 

2  .101,07  .UO./iOW  .000, ’  .2ti.i024  1.00  77  -.0340 

3  . 10  103  .044601  .0733  .250400  -1.5436  -4.8020 

4  .  19730  .040/97  -.0200  .032006  -2.3  1  44  -0.2392 


thetam  dthetam 


.0942 

.0964 

.0944 
.09  39 


- .0003 
.0022 
-  .0020 
-  .0005 


126 


r 

.  f.  n 

b 

.  v>  .  ii 

.11  .  o 

7, 

.499.::. 

. 09''.'  /;■ 

8 

.  50 j  .)/ 

.1-6.  40 

9 

.64215 

.  1  4.:'  Cl 

10 

.73'  3b 

.if;;:-.' l .: 

1  1 

.  79  ,b 

.03:. )  1  3 

12 

.  79b  19 

.020113  - 

13 

.  070  ’£' 

.  0  tl  1  7  1  2 

14 

1  3  L 

.122102  - 

19 

1.13619 

. 1 39091 

"  ’  -  .  ° n  !  > 


.  3  ■  3  .* : 

- 

:lj.  /  J 

.01  03!  ) 

i- 

.  2  0.3 4 

01  0  9 

-.2b’..  3  -  J 

•  <.  1  ae 

0.  .  9 

.291 391 

-1 

.  1  9  3  ! 

9  /  1 

- .2/ 1  -73 

- 

.  /  .  •  1  3 

UP  /0 

- . 244  320 

1 

•  tl  0 1 1  4 

0  17  9 

-  .  1  3  5  9  4  1 

■t 

.  0  3  3  9 

0  7  70 

.  26b'.  2U 

- 

.3123 

73l 

-.204300 

.1111 

0G1S 

.230391 

- 

.  9  4  3  o 

•  c 

• » i  •  r 

-  .  |M24 

.  f,.  '1 

. ')  j-  l 

-.0015 

5  .  7  02b. 

.  I.l  r; : .  4 

-.0017 

3  .  7  9  2  J 

.  9  3/  2 

-.0011 

9 . 9.  9  3 

.  i)0. -9 

-  .0013 

.9.11 

.  J  9 

.0010 

3.64 

.  •  '  q  t  ‘ 

-  .00 1 4 

A  .  9 9 6 

.0316 

- .0000 

4.329b 

.07  35 

-  ,0021 

3.9.03 

.0777 

-  .01'  l  8 

3.1437 

.0741 

-  .0037 

Impact"  1  errmom*  .  27939G772o-07  errvnl*  -  .  7  1  5255 7 3 7o-06 
t mpact*  2  errmom*  .  1  1  1  758709C-07  eri'vo)»  -  .  4  730  3 7  1  GGo-x/b 
impact*  3  errraeif.*  -.7450580000-08  orrvol*  .  33  0 2  7  .!C  9e -06 

impact.*  4  errmom*  -  .  35273U4 23c -07  orrvol*  .000p000U0e+00 

Impact*  5  errmom*  .  2 2 35 1  7 4  1  .,e-07  orrvol*  .  OtVv  .i0Oi'0e+00 

impact*  6  err  mom*  -.3725290300-00  orrvol*  .  OP  >900  .'i  3'o+00 

Impact*  7  or  rnoin*  -  .  2  3  20503  4  4  e-00  orrvol*  .  35 t  2  7  Op'jo-OO 

Impact*  0  errmom*  .  37252903;«v-X,'U  orrvcl*  . >.•>.>  i>  ‘U  )h 0o  +  00 

Impact*  9  or  mom*  -  .  7-1  503001.  ao-t  II  orrvol*  .  Oi  ;  •  :Ulnl’i)oUl3 

impact*  10  orri.icm*  —  .  1  1  1 /5O".'P337-07  orrvel*  .00.0>>DiM  000+00 
In. poet*  11  o.fi.iom*  .  9 3  l  32 29  / 5.. -09  crrvo-1*  -.1192092900-06 
Impact*  12  crrniuiii*  .4  912  721x30-07  orrvol*  .0000000000+00 
Impact-  13  errmom*  -  .  1  Of.7040  1  tie-03  orrvol*  -  .  23U4  1  U5/9e-06 
Impici'.  14  ori'i.iOni*  .  3  729290  J0<. -00  crrvol*  .  ilt'uiUij.1.11  00O+00 
impact*  IS  orrmom*  -  .  1  G763U0G3c--07  orrvol*  -  .  23U4 1 U579e-06 


1  mpact 

thotam 

t  imo 

dthetam 

dt  1  mn 

2 

.0964 

.  12  3  0 

-  .0035 

.1238 

S 

.09)5 

.3099 

-.004  9 

.2462 

8 

.0072 

.0177 

-  .004  3 

.  2470 

12 

.0310 

.  0  0  5  G 

- . 00 50 

.2479 

14 

.0777 

1.1137 

-  .0h'39 

.2481 

17 

.0727 

1.3017 

-.0050 

.2481 

20 

.0684 

1 . G09U 

- .0043 

.2400 

23 

.0040 

1 .8502 

-.0044 

.2404 

26 

.0594 

2  .  1000 

- .0046 

.2484 

28 

.0562 

2.3553 

- .0032 

.2407 

31 

.0512 

2 . 0035 

- .0050 

.2482 

33 

.048  2 

2.3525 

- .0030 

.2490 

36 

.0434 

3.100  1 

- .0049 

.2476 

39 

.0  395 

3 . 3439 

-.0039 

.2480 

43 

.0315 

3 . 1.!  4  5  9 

-.00133 

.4970 

45 

.0279 

4 .0959 

- .0037 

.2500 

47 

.0243 

4 . 3403 

- .0030 

.2505 

49 

.0206 

4.5973 

-  .0030 

.2509 

51 

.0170) 

4.0487 

- .0030 

.2514 

03 

.0133 

5  .  1  005 

- .0037 

.2518 

55 

.009  4 

5.3519 

-  .0033 

.2514 

57 

.  C056 

5.5972 

- .0039 

.  2453 

up  to  the  22 

peal  the  least  squares  fit  to  the  peaks  is  thetam»q+s*t  with  q* 

.  9574  1  tlOUCc-V  i  and  s*  -  .  1  C4  i  720540-01  with  wax  error*  .  2 72 37  1  1  8 1 6-02  at  peak* 
1  and  variance*  . 372083668e-04 


Impact  time  del  tat  thetas  x  tdotv  xdot  thetam  dthetam 


**“*■* 

“  *  . . 

. 

59 

5.80442 

.  199719 

.0009 

- .094315 

-.1293 

.0390 

.0051 

.0006 

60 

7.125271 

.24005-1 

-.0624 

.090963 

- .0981 

-.3624 

.004  5 

- .0006 

61 

7.62555 

.500273 

-  .003G 

- . 104315 

-.0690 

-.0549 

.0045 

.0000 

62 

7.66062 

.035070 

-.0043 

-.105942 

.0280 

.1148 

.6/045 

.  0000 

63 

7 . 670  69 

.019063 

- .0033 

—  .  1 03606 

.0/66 

.1933 

.0044 

.  0000 

64 

7.69195 

.012266 

-  .0622 

-.101259 

.0969 

.2209 

.004  4 

.0006/ 

65 

7.70070 

.000750 

-  .00 1 3 

- .099239 

.  1090 

.2451 

.0044 

.06)00 

66 

7 . 70803 

.000 1 25 

-.XU  U 

-.097233 

.1136 

.2519 

.0044 

.06100 

67 

8.33603 

.  62C007 

- .00 1 0 

.092121 

-.0063 

- . 4279 

.0030 

- .0006 

68 

8 .73245 

.  39561  6 

.0037 

-  .0:10246 

.0090 

.2501 

.63037 

-.0~02 

GO 

9 . 32574 

.595200 

- .0032 

.089229 

-.0311 

-  .  240  4 

.6'034 

- .0002 

70 

9.95603 

.  6  36'  2  d  3 

.0032 

-  .009240 

- .0269 

.0345 

.  00  3  4 

-  .0660  1 

71 14.95i 055 

.000020 

-  .  00  1  4 

.  0U832  7 

.  0000 

.0000 

.  0000 

-  .  066  3  4 

72 

.0>.’6’00**  **  *  *  *  * 

.0600 

.  000000 

.  001*0 

.0000 

.  66000 

.0)600 

73 

.00000 

.000000 

.06ii0 

.  0616,000 

.0000 

.  0000 

.0000 

.0000 

74 

.000  .10 

.06  pi  >.00 

.00  JW 

.001 000 

.0010 

,  0616*0 

.  X/00 0 

.0000 

75 

,00c  00 

.  00  0.0  00 

.06U0 

.000000 

.  0000 

.0000 

.006)0 

.0000 

76 

.  001  110 

.0006 j0 

.06' 0)1 

.000000 

.0060 

.000  0 

.  00610 

.0000 

77 

.01 M  (0 

.0h.-‘i  00 

.06  ;« 

.000000 

,  0000 

.0000 

.  001. 0 

.00/6  0 

70 

.01  U 

.00. 600 

.000:1 

.000000 

.00610 

.  000  0 

.  6661)161 

.  0)  *  6 1 0 

79 

.60:  v.iJ 

.  0»  r;'0U 

.0)  O0 

.060 U'U 

.  006-0 

.  0IO.I  0 

.  0000 

.0  >666  61 

80 

.8*  1  >)aj 

.  00  1-60 

.  -is  fil 

.006,  “.00 

.  6100  0 

.  6* )  6)  0 

.  JJ  »  0  66 

.0  J0  0 

01 

.00  0'J 

.06  0.-  •  0 

.  0-o-y 

.006000 

..06100 

•  0  60J  0 

.01100 

.0000 

■»  -» 

1. 

.00200 

.  0.'.  .  u< 

.  6  n  l  ■  0 

.  0  0  0  0  u  0 

.  0  0 1;  • 

.  06)00 

.  6.1000 

.00110 

83 

.00*.  60 

.  .06' 01.  ...  i, 

.6000 

.000000 

.0000 

.0000 

.  66  0  66  0 

.0000 

84 

.00000 

.060000 

.06U.10 

.006000 

.0000 

.0000 

.0000 

.0000 

Impact 

thetam 

1 1  me 

dthetam 

dt  1  mo 

66 

.6044  7 

.7729 

-  .0O5o, 

7  .  7  729 

* 

60 

.0037  0 

.  7364 

-  .006  7 

.9635 

up  to  the  2 

pc ik  the  least  squares  fit  to  the  peaks  Is  thetam=q+s*t  with  q* 
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